ここ で は , ミッ クス ト ・ シ グ ナ ル (アナ ログ ・ デ ィ ジ タル 混在 ) 
回 路 を , Verilog HDL シミ ュ レ ー タ を 使 つ て 検証 する 手法 を 
紹介 する . アナ ログ HDL な ど に 対応 する 高価 な シミ ュ レ ー タ を 
用 意 し な く て も , 低 コ スト の 環境 で 回 路 の ふる まい を 確認 で き 
る . サン プル は , トラ ン シ ー バ LSI な ど に 使わ れ て いる クロ ッ 
ク ・ デ ー タ ・ リ カバ リ 回 路 で ある . (編集 部 ) 


一 般 的 な ミッ クス ト ・ シ グ ナ ル LSI の 機能 検証 手法 を 議 
論 す る 場合 , 対象 と な る 回 路 に よっ て 最適 と 思わ れる 検証 
手段 が 異な り ま す . 例え ば , AD 変換 や D-A 変換 , PLL 
( phase-locked loop) 回 路 , PLL 回 路 を 含ん だ 機能 ブロ ッ 
タグ ( 今回 , 例 と し て 挙げ ば た クロ ッ ク ・ デ ー タ ・ リ カバ パリ 回 
路 や 周波 数 シン セ サ イ ザ な ど ), 電源 回 路 。 RF フロ ント エ 
ンド を 含ん だ 機能 ブロ ッ ク ( ミキ サ 回 路 を 含ん だ 機能 ブ 
ロッ ク な ど ), LCD ド ライ バ 回 路 の 検証 に お いて , 一 つの 
手法 だ け で 対処 する の は 困難 で す . た だ , いずれ の 分 野 も 
検証 時 間 の か か る アナ ログ 部 を どう 取り 扱う か が 一 つの ポ 
イン トド に だ なり ます 。 


デー タ 図 一 一 - 1 組 の 差 動 較 
ペア 較 1 

4 デー タ 較 本 誠 
クロ ッ ク 図 一 | 5) 較  1 


クロ ッ ク 情 報 を 図 1 個 の スト リー ム 図 受信 デー タ か ら 図 
含ん だ デー タ を 送信 凶 と し て 転送 較 クロ ッ ク を 再生 図 


図 1 CDR に よる デー タ の 送受 信 リ 


CDR に よる デー タ の 送受 信 で は , 受信 デー タ か ら ク ロッ ク を 再生 する . IC 外 


の クロ ッ ク ・ ス キュ ー の 問題 を 回 避 し , 高速 デー タ 転送 は 実現 で きる . 


ミッ クス ト ・ シ グ ナ ル , アナ ログ ・ デ ィ ジ タル 混在 , Verilog HDL, 


在原 栄一 


今回 は PLL を ベー ス に し た クロ ッ ク ・ デ ー タ ・ リ カバ パリ 
回 路 を 例 と し て , Verilog HDL を 用 いた 機能 検証 を 行い ま 
す . Verilog HDL を アナ ログ 機能 ブロ ッ ク の ビヘイビア ・ 
モデ ル と し て 採用 する 利点 と し て , MATLAB や 言語 の 
モデ ル よ り も 最終 的 な レイア ウト 情報 に 近い こと が 挙げ ら 
れ ま す . AD コン バー タ や D-A コン バー タ の よう な 離散 時 
間 的 な アナ ログ 機能 ブロ ッ ク の ビヘイビア ・ モ デル の 記述 
に Verilog HDL を 用 いた ケー ス は 珍し く あ り ま せん で し 
た . し か し , PLL の よう な 連続 時 間 的 な アナ ログ 機能 
ロッ ク の ビヘイビア ・ モ デル を Verilog HDL で 記述 する 例 
は , それ ほど 多く あり ませ ん で し た. 

以降 で は , 対象 と な る クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 
を 紹介 し , 検証 方 法 や ビヘイビア ・ モ デル の 内 容 , シミ ュ 
レー ショ ン 結果 に つい て 解説 し ます . 


4 クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 


現在 の 高 性 能 通信 シス テム で は , デー タ 重り 1Gbps 以 
上 の 送信 速度 が 要求 され る ケー ス が し ば あり ます . 


人 @ 高速 デー タ 伝 送 で は クロ ッ ク ・ デ ー タ ・ リ カバ リ か 必須 

一 般 的 な シン グル エン ド 標準 1/O 規 格 で は , 250MHz 近 
辺 の 周波 数 で ノイ ズ 制 約 に 直面 する と 言わ れ て いま す . こ 
れ に 対し て 差 動 標 準 1/O 規 格 で は , コモ ン ・ モ ー ド の 除去 
に よっ て , より 少な い ピ ン 数 で 高速 な デー タ 送信 を 可能 に 
し まし た . し か し , 差 動 標準 {/O 規 格 で も , 1Gbps 近 くに 
な る と , クロ ッ ク ・ ス キュ ー が 問題 と な っ て きま す . 


回 jmBL。 多 回 = デー タタ ュ 萌 交 0 』 


CDR, VeritakWin, M 系 列 , 疑似 乱数 。 ペ ルン ゲ ク ッ タ 法 
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クロ ッ ク ・ デ ー タ ・ リ カバ リ ( CDR : clock data recov- 
ery) は , 受信 側 に お いて デー タ 列 か ら ク ロッ ク を 再生 する 
方 式 で あり , 受信 側 IC 外 の クロ ッ ク ・ ス キュ ー の 問題 を 
避 し , デー タ と クロ ッ ク の 同期 を 容易 に し まず 図 1). そ 
の 反面 , CDR で は デー タ が 基本 的 に NRZ non return 
zero) 形式 で ある た め , デー タ 列 か ら VCO voltage 
controlled oscilator : 電圧 制御 発振 器 ) の 適正 な 発振 中 心 
周波 数 を 求め る の に リフ ァ レ ンス 周波 数 が 必要 と な り ま す . 
も し く は , DLI( delay locked loop) な ど を 使用 し た より 複 
雑 な 回 路 を 付加 し な けれ ば な り ま せん . 

な お , CDR で は , デー タ お よび 内 部 生成 クロ ッ ク の ジッ 
タ 量 例え ば , 外来 ノイ ズ や , ルー プ の 持つ 周期 性 や デバ 
イス ・ ノ イズ な ど に 起因 する ノイ ズ ) や ロッ ク す る まで の 
時 間 な ど が , 性 能 を 決め る 指標 に な り ま す . また , ユー ザ 
側 か ら 見 る と , リフ ァ レ ンス ・ ク ロッ ク が 不要 に な る , 広 
い 範 囲 の デー タ ・ レ ー ト ( クロ ッ ク 周 波数 ) に 対応 で きる , 
ロー パス ・ フィルタ ( LPF : low pass filter ) 用 の 外 付け 
キャ パシ タ が 不要 に な る , VCO を 内蔵 する , と いっ た 特 
徴 が 望ま れ ま ず 現在 , 半導体 メー カカ から 提供 され て いる 
多く の カス タム 品 で , 実現 され つつ ある ). 


回 


@ 10Gbps に 対応 する サン プル 回 路 

CDR は , 現在 , 比較 的 多く の 半導体 メー カ か ら カ スタ ム 
品 , も し く は 機能 ブロ ッ ク と し て 提供 され て いま す . ここ 
で は , CDR の アー キテ クチ ャ と 回 路 に つい て , 数 年 前 に 
IEEE Solid State Circuits 誌 に 掲載 され た 論文 を 例 に 説明 


クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 


ディ ジタル ・ シ ミュ レー タ で 検証 する 


し まず うめ . 仕様 は 10Gbp& クロ ッ ク 周 波数 は 5GHz) で す . 
使用 プロ セス は 018m プロ セス の CMOS 技 術 で , VCO を 
内 蔵 し て いま 実際 の 製品 と し て 考え た 場合 , レイ アウ 
ト 設計 へ の 負担 も 含め て , と て も 難し い 仕 様 の よう に 思わ 
れる ). この CDR の 例 を , 今回 の 機能 検証 の 対象 と し ます . 

図 2 が CDR の アー キテ クチ ャ で す . VCO の 周波 数 を 高 
くし た く な い の で , 10Gbps の データ ・ レー ト を 実現 する 
の に 5GHz の クロ ッ ク の 半周 期 に 対し て 1 個 の デー タ を 割 
り 当て て いま す . また , ロッ ク 可能 な 周波 数 範囲 を 広げ る 
目的 ,」 お よび ロッ ク す る まで の 時 間 を 短縮 する 目的 で , 位 
相 検出 ・ 比 較 回 路 ここ で は PD と 呼ぶ ) だ け で な く , 周波 


受信 デー タ 

( 10Gbps) 較 
ハー フレ ー ト 較 ハー フレ ー ト 較 Sa 
FD PD ( 10Gbps) 


電圧 - 電 流 図 
コン バー タ 罰 


再生 クロ ッ パ 5GHz) 


図 2 CDR の アー キテ クチ ャ 例 
5GHz ク ロッ ク で 10Gbps の スル ー プ ッ ト を 得 ら れる アー キテ クチ ャ と な っ て 


1、 1 


1 


ゆみ 


B 


Em: 


Hg] 


2 


( a) 0", 45", 90「, 


135* の 位相 差 を 持つ 正弦 波 電圧 生成 回 路 較 


Vg 


ーーーO 7 Oーーーー ゅ \C 
- | オーー 
の | 間 Muze JI 

図 3 1 ゆめ 
VCO チ ェ ー ン 回 路 /ss 
4 個 の VCO を 縦列 接続 する こと で , 
0t deg), 45" , 90", 135" の 位相 回 ( 
差 を 持つ 正弦 波 電圧 を 生成 し て いる . ( c) VCO 回 路 図 ( b)( a) の 各 ノ ー ド に 対す る 電圧 【 
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数 検出 ・ 比 較 回 路 FD と 呼ぶ ) を 設け て いま す . さら に , 
VCO で 生成 する 4 種類 の 位相 攻 0', 45 , 90", 135") 
を 持つ クロ ッ ク を 使用 し , PD や EFD の 内 部 で , デー タ の 持 
つ ク ロッ ク 情 報 と VCO で 生成 し た クロ ッ ク と の 位相 進み 
や 位相 遅れ , 周波 数 の 大 小 を 判断 し ます . 

3 が , 4 種類 の 位相 差 を 持つ クロ ッ ク を 生成 する VCO 
チェ ー ン 回 路 で す . 図 a) は , 図 3 c) の よう に 4 個 の 7C 
タン ク 型 VCO を チェ ー ン 接続 し た 回 路 で す . 四 つ の ノー 
ド C A, B, C, D) で 45" ずつ の 位相 差 の 正弦 波 電 圧 を 生成 
し まず 図 3 b)). 

同じ よう に , 図 4 は PD 回 路 お よび FD 回 路 で す . PD で 
は , 受信 デー タ に 対し て 内 部 生成 クロ ッ ク の 位相 が 進ん で 


( a) PD 回 路 較 


図 4 PD 回 路 と FD 回路 
4 種類 の 位相 差 の クロ ッ ク を 用 いて , 位相 差 や 周波 数 差 を 検出 する . 


内 部 生成 クロ ッ ク の 位相 が 遅れ る 場合 較 


S+ 図 

「 区 図 i 1 
Din ! 1 ! ! 1 
生ま 
内 部 生成 クロ ッ ク の 位相 が 進む 場合 図 

ゞ - 図 

Dn 9 和 il 遇 間 
ce ルキ キル ルト TT 


a) 受信 デー タ と 内 部 生成 クロ ッ ク と の 位相 関係 の 検出 図 


図 5 受信 デー タ と 内 部 生成 クロ ッ ク の 位相 。 お よび 周波 数 の 関係 の 検出 


Vspi は 2 値 を , Vp は 3 値 を 表現 する . 
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クロ ッ ク 周 波数 が 小さ い 場 合 


いる か 遅れ て いる か の 判定 を 行い ます . FD で は , 2 個 の 
PD に 45” ずつ 位相 差 の 異な る (VCO に よる ) 内 部 生成 ク 
ロッ ク を 入力 し , 周波 数 の 大 小 の 判定 を 行っ て いま す . こ 
れ ら の 判定 の よう す に つい て は , 図 5 を 参照 し て くだ さい . 
図 & a) で 受信 デー タ D。 の 論理 状態 に 変化 が な い 場合 は , 
判定 を 保留 し て いま す . また , 図 呈 b) で は , 周波 数 の 大 小 
判定 を Vp と Vpps の 位相 差 で 決め て いま り ぁ の 符号 ). 

図 6 は チャ ー ジ ・ ポ ンプ 回 路 ここ で は CP と 呼ぶ ) で , 
図 人 4 b) の ゆり j と っ を 入力 と し て 切 ー を 出力 し ま 
す . な お , 出力 pn- V ぁ に つい て は , 図 & b) も 参照 く 
だ さい 、. 


図 7 は ロー パス ・ フ ィ ル タ で す . 図 6 の チャ ー ジ ・ ポ ン 


( b) FD 回 路 較 


クロ ッ ク 周 波数 が 大 きい 場合 凶 


Wepn- Vp し 1 [| 1]. Mnp-Wp し レー L 


MPp1 ーー トーー ト ーー オ - WPp1 ーー トーー ト ーー 
Vpz ー ps 一 テ 


Up ニーーー ト ーー 1 生ま 


( b) 受信 デー タ と 内 部 生成 クロ ッ ク の 周波 数 の 関係 の 検出 図 


いい 


プ 回 路 出 力 で ある Vp - Vp ゅ が この 回 路 の 入力 。 に 接続 さ 
れ ま す . さら に , この は VCO の 制御 入力 と な り , 
3 c) の し, と 接続 され ます . この フィ ル タ に よっ て PLI( 二 
つの 極 と 一 つの 零 点 を 持つ) を 構成 し ます . な お , Cs は ク 
ロッ ク ・ フ ィ ー ル ド ・ ス ルー に よる 電圧 リ プル を 抑制 する 
た め の 容 量 で , 図 7 の 構造 は と て も 一 般 的 で す . 


『: 機能 検証 の 方 法 


現在 , 比較 的 大 規模 な モデ ル , また は 機能 検証 時 に 多く 
の イベ ント を 要する モデ ル の ミッ クス ト ・ シ グ ナ ル ・ シ ミュ 
レー ショ ン に は , CO 言語 や MATLAB, Verilog-A, 
Verilog-AMS な ど が 用 いら れ て いま す . シス テム 検証 か ら 
LSI の 作り 込み まで の どこ に 焦点 を 当て る か に よっ て , こ 
れ ら の 言語 や ツー ル が 使い 分 けら れ て いま す , 


@ ディ ジタル ・ シ ミュ レー タ で アナ ログ 要素 の 動作 を 模擬 
今回 は , Verilog HDI( Verilog-A と 対比 され て , 
Veriog-D と も 呼ぶ ) 言語 に よっ て , 上 記 で 説明 し た CDR 
路 の 機能 検証 ミッ クス ト ・ シ グ ナ ル ・ シ ミュ レー ショ ン ) 
を 行い ます . ミッ クス ト ・ シグナル 回 路 を ディ ジタル 回 路 
用 の Verilog HDL に よっ て 検証 する 意義 は , 以下 の 3 点 で 
ず 図 8). な お , Verilog HDL を コン パイ ル ・ 実行 する シ 
ミュ レー タ と し て , ここ で は Windows XP 上 で 動作 する , 
菅原 シス テム ズ の VeritakWin を 使用 し まし だ 下 掲 の コラ 

人 VeritakWin と は ? 」 を 参照 ). 

e 合成 対象 と な る ディ ジタル 回 路 は RTI( register transfer 
level) 記述 を 使っ て 機能 検証 を 行い た い た だ し , 今回 の 
例 で は , 合成 対象 と な る ディ ジタル 回 路 は 多く な い ). 

e アナ ログ 回 路 に お いて も , LSI の 作り 込み に 近い モデ ル 
を 使っ て 機能 検証 を 行い た い . つま り , レイ アウ ト ( で 
き 上 が っ た LS1I) と 同じ 端子 情報 を 備え た モデ ル で , 機能 
検証 を 行い た V い ( Verilog HDL に よる モデ リン グ で は 
アナ ログ 端子 を 表現 する た め に , 要求 され る 精度 に よっ 
て バス 幅 が 決ま る . 例え ば ANALOG_IN と いう 端子 は , 
Verilog HDL で ANALOG IN[15: 0] と 表現 する ). 

e 用 意 し や すい 多く の ユー ザ を 有する ) 設計 環 境 で , ミッ 
クス ト ・ シ グ ナ ル LSI の 機能 検証 を 短 時 間 で 行い , 全体 
性 能 の 向上 を 妨げ る 要因 アー キテ クチ ャ や 回 路 ) を 特定 
し た い イベ ント 数 を 最小 限 に と ど め , 処理 時 間 の か か 
る ニュ ー ト ン ・ ラ プ ソ ン 計算 な どの 数 値 演算 を で きる だ 


回 


クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 
ディ ジタル ・ シ ミュ レー タ で 検証 する 


ニー ニー 
PDi か ら 較 - | 央 過 FD か ら 
M1 M2 Ms M4 
/2 
図 6 チャ ー ジ ・ ポ ンプ 回 路 
電圧 - 電 流 変換 回 路 を 兼ね て いる . 
訪 O や 
x1 
還 pp 上 cz 
図 7 0 
ロー パス ・ フ ィ ル タ 〒 
PLL を 構成 し て いる . x1 と x2 は 状態 方 
程 式 で 表現 し た 場合 の 状態 変数 と な る 
ノー ド を 示し て いる . V 


け 抑 える こと に より , SPICE や Verilog-A と 比べ て 検証 
時 間 を 短縮 する ). 


@ ジッ タ 量 の ピー ク 値 と RMS 値 で 評価 
ここ で は , も っ と も 代表 的 な 特性 値 で ある ジッ タ 量 の ピ 
ー ク 値 と RMS root mean square) 値 を 評価 基準 と し ます . 
た だ し , ロッ パ ( 収束 ) す る まで の 時 間 を 1.5s 以 内 と し ま 
ポ つま り , ジッ タ 量 の 測定 を 開始 後 1.54s 後 か ら 始 め る ). 
今 , SONETK synchronous optical network : 基幹 伝送 
網 の 多重 化 通信 体系 ) の 規格 に 沿っ て , ピー ク 値 を 01UI 


《VeritakWin と は ?》 


VeritakWin は, 菅原 シス テム ズ http://japanese.sugawara- 
syStems.com/ ) の Verilog HDL シミュレー タ で , 編集 , デバ ッ 
グ , 波形 観測 の 機能 を 有 し て いま す . Windows XP と Windows 
2000 の 上 で 動作 し ます . な お , 筆者 が 使用 し た バー ジョ ン は 


1.71A で す . 

本 シミ ュ レ ー タ の ユニ ー ク な 特徴 と し て , 時 間 連 続 系 で 常 微 
分 方 程 式 を 解く た め の ラ イブ ラリ を 用 意 し て いま す . SPICE な 
どの アナ ログ ・ シ ミュ レー タ の 過渡 解析 に 相当 する シミ ュ レ ー 
ショ ン を 行え ます . 
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二 に 計 


軒 出 峰 潤 


回 路 イ メー ジ を 想定 し , 劣化 還 
要因 素子 ば ら つ き や 非 線形 
要素 な ど ) の 列 拳 を 行う 凶 


アー キテ クチ ヤ ( 方 式 と アル ゴリ ズム ) を 想定 図 


全体 を ファ ンク ショ ン ・ ブ ロッ ク に て 実現 図 


最終 的 な 端子 情報 を 較 
…③ 図 \ 反映 し た も の で は な く 図 
て も よい レイ アウ ト 罰 
と の 一 致 は 問わ な い ) 図 


…@ 図 1 
/ 


Verilog HDL に よる 検証 は , 較 


回 路 実 現 を 行い , 劣化 要因 図 
( 素子 ば ら つ き や 非 線形 要素 了 
な ど ) の 列挙 を 行う 凶 


…⑤ 較 上 下 の フ ロー で 下 側 に 近い 図 

ー 場所 に 位置 する . 簡便 な 設計 図 
環境 で , レイ アウ ト に 近い 図 
端子 情報 を 持っ た ミッ クス ト ・ 匠 
シグ ナル ・ シ ミュ レー ショ ン 罰 


…⑥ 図 を 比較 的 短 時 間 で 行う こと が 鐘 
で きる 図 
…⑦ ーー ン バ 
最終 的 な 端子 情報 を 較 


( 回 路 に よる 劣化 要因 を 考 


シミ ュ レ ーション の 実行 Veriog-AMS な ど ) 図 …⑨ 


回 路 ブ ロッ ク に 対応 する ビヘイビア ・ モ デル の 作成 
慮 し た モデ ル が 必要 ) 囚 


…⑥⑧ 図 \ 正確 に に 反映 し た も の 較 
で ある こ @( レイ アウ ト 図 
と は 完全 に 一 致 ) 図 


設計 の 終了 図 


図 8 ミッ クス ト ・ シ グ ナ ル LSI の 設計 ・ 機能 検 証 フ ロー 例 と Verilog HDL を 使用 する 意 
アナ ログ ・ ビ ヘイ ビア ・ モ デル を 用 いた 設計 フロ ー の 例 を 示す . SPICE 系 の アナ ログ ・ シ ミュ レー タ で は 多く の 解析 時 間 を 要する た め , ビヘイビア ・ モ デル の 尊 


入 を 行い , 段階 を 踏ん だ 階層 設計 ・ 検 証 が 現実 的 な 解 と な る . 


( unit intervals) 以内 , RMS 値 を 001UI 以内 と し ます . 今 
回 , UI = 100ps と する の で , ピー ク 値 は 10ps が , RMS 値 
は 1ps が 目標 と な り ま す . な お , ピー ク 値 と は 偏差 の 最大 
値 と 最小 値 の 差 , RMS 値 と は 偏差 の 2 乗 和 の 平方 根 を 求め 
た も の で す . 


デー タ 列 を 生成 する 較 
リフ ァ レ ンス ・ ク ロッ ク 


図 9 ジッ タ 評 価 の た め の ス キュ -( 時 間 差 ) の 測定 
生成 クロ ッ ク ・ エ ッ ジ か ら リ ファ レン ス ・ ク ロッ ク ・ エ ッ ジ まで の 時 間 差 の 
変動 を ジッ タ 量 と する . 
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また , ジッ タ 量 の 規定 で す が , 入力 デー タ 列 を も と に 3 
成 し た クロ ッ ク ・ エッ ジ か ら , デー タ 列 を 生成 する リフ ァ 
レン ス ・ ク ロッ ク ま で の 時 間 差 の 変動 を 測定 し まず 図 9). 


『: モジ ュー ル の 記 可 


ここ で は シミ ュ レ ーション を 行う テス ト ベン チ , CDR 回 
路 , 入力 デー タ 生成 回 路 , 測定 回 路 の 順 で , それ ぞ れ の 記 
述 を 紹介 し ます . 


EE 


⑯ テス ト ベ ンチ : デー タ 生 成 と ジッ タ 測 定 の 機能 を 用 意 

10 が テス ト ベン tb_topvy) の ブロ ッ ク 図 で す . 入力 
デー タ 生成 回 路 dat_genv) で 生成 され た デー タ 列 Data) 
は , CDR 回 路 topvY) へ と 入力 され , 内 部 で 生成 され た ク 


クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 『 
ディ ジタル ・ シ ミュ レー タ で 検証 する 


《M 系 列 に よる 疑似 乱数 の 生成 》 


M 系 列 信 哨 maximum length sequence) は 2 値 “ L" ま た ば H”) 
の 周期 信号 で す . 具体 的 な 回 路 は 図 N-1 の と お り で , ヵ 個 の DFR D 
フリ ッ プ フロ ッ プ ) と 最低 1 個 の XOR 排他 的 論理 和 ) ゲ ー ト で 最大 
2 1 長 の 乱数 を 生成 で きま す . また 表 A-1 は , ヵ の 具体 値 と , 最大 
2%- 1 長 の 乱数 を 得る た め の XOR 入力 の タッ プ 位 置 を 示し て いま す . 
2%ー 1 長 の 乱数 が 繰り 返さ れる こと か ら , 疑似 乱数 列 と 呼ば れ て い 
まず め . 


図 A-1 シフ ト ・ レ ジス タ 罰 
M 系 列 信号 の 
生成 回 路 

思 個 の DFF と 最 
低 1 個 の XOR 
で 最大 2 の - 1 長 
の M 系 列 信号 


表 A-1 M 系 列 信号 に お ける 最大 系 列 長 と XOR 
入力 タッ プ 位 置 
シフ ト ・ レ ジス タ 長 M 系 列 長 
( 記 ) 語り 
15 32767 
131071 
262143 
1048575 
2097151 
4194303 
M 系 列 信号 図 8388607 
01101001… コラ 
268435455 
536870911 
2147483647 


( 乱数 ) を 生成 で 
きる 


8589934591 


リフ ァ レ ンス ・ 


リフ ァ レ ンス ・ 図 EYE 


クロ ッ ク 図 


, クロ ッ ク 図 


CP 処 理 , ロー パス ・ フ ィ ル タ 処 理 較 | vco 処 理 図 
( cp_lpfy) 較 ( vco_chnv) 良 一 =-| ( msjv) 


デー タ 生 成 図 | ! | : ーーーーーーー コ 
23 ビ ッ ト の M 系 列 隊 ! ! で 人 
( qat_genv) 較 4 ! | gtl_dly.v) 凶 


テッ 


図 10 CDR 回 路 top.v) の テス ト ベン チ 


ー L  , ピー ク /RMSE 
デー タ 再 生成 還 ! ジッ タ 図 
( re_time.v) 凶 


toDp.v 


再生 成 さ れ た デー タ 図 


テス ト ベン チ を 示す . 大 きく トッ プ 回 路 CDR 回 路 ), 入力 デー タ 生成 回 路 。 ジッ タ 測 定 回 路 か ら 構成 され る . 


ロッ ク で 読み 込ま れ , 出 ガ Re-timed DataJ) さ れ ま す . ま 
た , 測定 回 路 msjty) に よっ て , VCO 回 路 vco_chnY) 
で 生成 し た クロ ッ ク と , 入力 デー タ 列 を 生成 し た り ファ レ 
ンス ・ ク ロッ ク の 時 間 差 が 測定 され , ジッ タ 量 が 求め られ 
ます , 

リス ト 1 が テス ト ベン tb_topy) の ソー ス ・ コ ー ド で 
す . 各回 路 で 使用 する パラ メー タ は , parameter 文 に 
よっ て ここ で 設定 し , qefparam 文 で 下 の 階 層 モ ジュ ー ル 
へ と 引き 渡し ます . パラ メー タ の 意味 に つい て は , 該当 す 
る モジ ュー ル の と ころ で 説明 し ます . また , 時 間 単 位 と 精 
度 は と も に 1ps と し ます . この 1ps は , ジッ タ 量 RMS 値 ) 
の 目標 に 対し て 余裕 を 持っ て は いま せん が , 測定 は 可能 で 
す . な お , コー ド の 簡単 な 説明 を リス ト 1 の 中 に コメ ント 
文 と し て 記し て あり ます . 


@ CDR 回 路 : 階層 化 さ れ た 五 つ の ブロ ッ ク か ら な る 

再度 , 図 10 を 参照 し て くだ さい . topyY が CDR 回 路 の 
トッ プ 階 層 top) です. この top モジ ュー ル と , top を 構成 
する モジ ュー ル module) に つい て , 以下 に 説明 し ます . 
1) top モ ジュ ー ル top.v) 

top モ ジュ ー ル は CDR 回 路 の トッ プ で , 図 2 に 対応 し て 
いま す . リス ト 2 p.121) が ソー ス ・ コ ー ド で す . テス ト ベ 
ンチ 上 で 指定 され た パラ メー タ 値 を parameter 文 で 受け 
取り , defFparam 文 に よっ て 下 の 階 層 モ ジュ ー ル へ と 引き 
渡し ます 。 パ ラメ ー タ の 意味 に つい て は 。 該当 する モ 
ジュ ー ル の と ころ で 説明 し ます . 
2) vco_chn モ ジュ ー ル ( vco_chn.v) 

vco_chn モ ジュ ー ル は VCO 回 路 VCO チ ェ ー ン 回 路 ) 
で , 図 3 a) を 実現 し た も の で す . リス ト 3 pp.122-123) が 
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リス ト 1 テス ト ベ バン tb_top.v) 


テス ト ベン チ の パラ メー タ 値 の うち , ディ ジタル 選 延 代 TD_PD1 と TD_FD) を 2 倍 に し て シミ ュ レ ーション し た 結果 . 図 15 に 対し て , ロッ ク 状 態 に お ける 大 き 


な リプ ル を 確認 で きる . 


//CDR 回 路 用 テス ト ベン チ 
//Revrsed on 01:29 05/11/23 


~timesca1e 1pg / 1pg 


~1nc1ude "Eop .Y リ 
~inc1ude "dat gen.Y リ 
~nc1ude "mg]t .Y リ 


modu1e 上 b op: 
reg エ SE 上 , ref co1k: 
Fed Str ヒ ms] 上 , okgen : 
re m8]t clk: 
re m8]t cl1k2: 
wire ou p, ou n, okgen45 p, ckgen45 n, da p, dat nz 
Wire YOd, VS8: 
wire [23:0] o1f ou: 


1nteger Ss refo1k: 
integer mcdz: 

integer tmp 11,EmD 12: 
rea1 て tmp 1,EmDp r2: 


// ト ッ プ 回 路 (top) と 測定 回 路 (ms]t) 用 パラ メー タ 
parameter TM PRCS=0.001 , / / [nsec] , 時 間 単 位 


// ト ッ プ 回 路 (top ) 用 パラ メー タ 
parameter K RTOT=524288.0: 
parameter TD PD1=0.3: // [ns] , 遅延 時 間 , pp 


parameter TRF PD1=0.1: // [ns] , 立ち 上 が り ・ 立 ち 下 が り 時間, pp 


// (TD+TRF) /K DrV 値 を 処理 サイ クル 
(cyc1e) と する . ED 
parameter TD FD=0.5: ng] , 遅延 時 間 , Fp 


parameter K DITV PD1=20: 


parameter TRF FD=0.1: ns] , 立ち 上 が り ・ 立 ち 下 が り 時 間 , Fp 


(TD+TRF) /K DIV 値 を 処理 サイ クル 
(cyc1le) と する . FD 


parameter K DTV _FD=20: 


parameter CUR CP PD1=125.0: 


//rea1 変数 か ら integer 変数 へ 変換 係数 


parameter RTO_B A PD1=1.0: 


// チ ャ ー ジ ・ ポ ンプ の 吸い 込み 電流 と 吐き 出 


電流 の 比率 (Un-match の 影響 度 ), PD 


uA ・ ポ ンプ の 電流 値 , pp 


Parameter CUR CP FD=125.0: // [uA], チャ ー ジ ・ ポ ンプ の 電流 値 , Ep 


parameter RTO B A FD=1.0, // チ ャ ー ジ ・ ポ ンプ の 吸い 込み 電流 と 吐き 
電流 の 比率 (On-match の 影響 度 ), rp 
parameter OFS CUR CP FD=0.0: // [uA], 中 間 レ ベル 入力 の 
チャ ー ジ ・ ポ ンプ の 電流 値 , Fp 
// [pF], rpPFP の キャ パシ タ 値 (フィ ル タ 定数 ) 
// [pF], チャ ー ジ ・ ポ ンプ 切り 替り 時 の 


parame セ ter CEP=50.0: 
parameter C2=2.0: 


フィ ー ド スル ーー 緩和 用 キャ パシ タ 値 . rp で 使 
// [Ohm], ロー パス ・ フ ィ ル タ の 抵抗 値 
(フィ ル タ 定数 ) 
parameter TS RK LPF=0.001: // [ns], ルン ゲ ク ッ タ 処理 の 周期 , rp 
parameter NN STP TNT LPF=0.064: // [ns], ルン ゲ ク ッ タ 計算 の ステ ッ プ 
初期 値 (rpr 出力 の 更新 周期 の 初期 値 ) 
// 許 容 相対 エラ ー 項 (ルン ゲ ク ッ タ 計算 ) 
// 許 容 絶対 エラ ー 項 (ルン ゲ ク ッ タ 計算) 


parameter RP=30.0: 


parameter RELERR TLPF=0.001: 
parameter ABSERR TLPF=0.00001: 
parameter CYCLE VCO=1 
parameter THETA TNITT=0.0: // [rad], VCO の 初期 位相 


// 測 定 回 路 (msjt) 用 パラ メー タ 
parameter CYCLE _MJ=1: 


// 測 定 を 行う 際 の スキ ッ プ 数 . 1/CycLE MT が 
測定 頻度 を 意味 する 

parameter T AVRG=300: // [ns], strt msjt が エ r->H 後 の ジッ タキ 

を 求め る た め の 時 間 を 指す 


出 


の 


//Vco の 出力 更新 時 間 (=TM PRCS*CYCLE VCO) 


F 均 値 


parameter T AVMJ=500: // [ns], ピー ク ・ ジ ッ タ 値 や RMs ジ ッ タ 値 を 


求め る た め の 時 間 を 指す 


parameter TS CLK=0.2: // [ns], リフ ァ レ ンス ・ ク ロッ ク の 周 


defFparam 12.TM PRCS=TM PRCS: 
defFparam 12.K RTOT=K RTOT : 

defFparam T2.TD PD1=TD PD1: 

defFparam T2.TRE PD1=TRE PD1: 
defEparam T2.K DIV PD1=K DTV PD1: 
defparam 12.TD EFD=TD_FD: 

defFparam 12.TRE FD=TRE _FD: 

defEparam 12.K DTV FD=K DIV FD: 
defEparam 12.N STP TNT LPF=N STP TNT LPF: 
defparam 12 .RELERR LPF=RELERR TLPF: 
defFparam 12 .ABSERR LPF=ABSERR TLPF: 
defFparam 12.TS RK LPF=TS RK LPF: 
defEparam 12.CUR CP PD1=CUR CP PD1 
defFparam 12.RTO _B A PD1=RTO B A PD1: 
defEparam T2.CUR CP FD=CUR CP FD: 
defFparam 12.RTO _B A FD=RTO B A FD: 
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deFparam T2.OFS CUR CP FD=OFS CUR CP FD: 
defFparam T2 .CP=CB: 

defFparam 2 .C2=C2: 

defparam 12.RP=RP: 

defparam 1T2.CYCLE VCO=CYCLE VCO: 
defparam T2 .THETA TNTT=THETA TNTT: 


defFparam T3.TM PRCS=TM PRCS: 
defFparam T3.CYCLE MJ=CYCLE MJ: 
defFparam T3 .TT AVRG=T AVRG: 
defFparam T3 .T AVMJ=T AVMJ : 
defFparam T3.TS CLK=TS CLK: 


//--- ネ ッ ト リ スト 

a881gn dd=1!D1 : 

a8819n Y88=1!b0: 

ag81gn ckgen45= (ckgen45 p==ckgen45 n)? 1'b0 
ckgen45 D: 

dat gen +T1 (rgt,ref oc1k,dat Dp,dat n): 

top 12 (rst, da _p, da n, ou p, ou n, 

ckgen45 p,okgen45 n, o1f ou 上 , vdd, YS8) : 
mg 上 3 (strt ms]t, ref co1k, okgen45) : 


initia1 begin 
tmp_r1=TS CLK/TM PRCS: 
tg_refo1k=tmD エ 1: 
// 入 力 デー タ 生成 用 リフ ァ レ ンス 
: 放 回 み タ 


mocd=$Fopen ( "rs1t . dat") : 
// 結 果 書 き 込み 用 ファ イル の 指定 
// リ セッ ト 有効 


58 ヒ =1「D1 : 
8 モエ ms] ヒ =1'D0: 


// 測 定 無効 


#100 rst=1'b0:  // リ セッ ト 解 除 
1200000 st て ms] ヒ =1'D1 : 
// 測 定 有効 
#810000 begin 
SEdigsp1ay (mcd, ッ タ 量 -- - 平 値 : sd",13.rs1t av) 
$Edisp1ay (mcd, "ジッ タ 量 -- - 正 側 ピー ク 値 : *q、 
負 側 ピー ク 値 : qr 
, エ 3 .m] pk max, 13 .m] pk min) : 
tmp エ 1=T3 .m] rms2 1: 
tmp_r2=T3 .hit trmsz 
$Edisp1ay (mcd, "ジッ タ 量 ---RMS の 2 乗 値 : sgr, 
(tmp_r1/tmp_r2) ) : 
$Edisp1ay (mcd, "VCO 生成 クロ ッ ク 半 周期 
ーーー 最 小 値 : *d, 
最大 値 : qr, エ 3.tsh min gen c1k, 
-tsh maxx gen Co1k) : 
$Edisp1ay (mcd, "ルン ゲ ク ッ タ 計算 で 許容 誤差 を 
満足 し な か っ た 回 数 : 
き を d, エ 2 .T2 .hit min h) : 


$Eoc1ose (mcd) : 
$Einigh: 
end 


end 


a1ways begin 
// 入 力 デー タ 生成 用 リ ファ レン ス ・ ク ロッ ク の 生 
ref co1k=1'b0: 拉 (tg refFo1k/2) : 
ref co1k=1'b1: 韻 (tg refFo1k/2) : 
end 


a1ways begin 
//Vco チ ェ ー ン 回 路 制 御 電圧 値 を 書き 込む た め の 
クロ ッ ク 生 成 


ms]t cl1k=1'b0: #500: 
ms]t cl1k=1'b1: #500: 
end 


a1way8 @ (ms]t cl1k) begin 
//Vco チ ェ ー ン 回 路 入 力 (制御 ) 電圧 値 の 書き 込み 
$Edisp1ay (mcd, "%d きり , き tme, T2.T3 .Y1n エ 1) : 
//1 行 目 は 時 間 , 2 行 目 は vco 入 力 


end 


endmodu1e 


ソー ス ・ コ ー ド で す . 機能 は リス ト 3 の 中 に も 簡単 な 説明 
を 記し て いま す が , 入力 電圧 を 図 1T p.123) の 関係 に 基 づ 
いで 4 次 多項式 近似 に より ) 45" ずつ 位相 の 異な る 四 つ の 
クロ ッ ク へ と 変換 し ます . また , パラ メー タ の TM PRCS 
は シミ ュ レ ーション 時 間 単 位 ~timescale で 与え た の と 
同じ 値 ), CYCLE VCO は スキ ッ プ 数 , CYCLE VCO* 
TM PRC8 値 は 出力 更新 期間 と な り ま す . それ か ら , 
THETA TNTT は 発振 出力 の 初期 位相 , K RTOT は 内 部 変数 
( 実数 ) と 信 哨 整数 表現 ) の 引 係数 ) で す . な お , 本 来 , 


リス ト 2 top モ ジュ ー ル top.v) 


//CDR 回 路 ( ト ッ プ 回 路 ) 
//Revised on 00:24 05/10/08 


~1no1ude "pFd top .Y リ 
~1no1ude "cp 1pFf .Y リ 
~1no1ude "voo chn. マ リ 
~1no1ude "re 上 1me . マ リ 


modu1e op (rst, dQdatp, datn, outp, outn, ckgenp, Ckgenn, 
C1fE out, vdd, YS8) : 
1nput rg, datp, datn: 
1nput YOd, YS8: 
Output outD, Outn: 
output ckgenp, ckgenn : / / 生成 クロ ッ ク で , デー タ を 生成 する 
リフ ァ レ ンス ・ ク ロッ ク と は 90” 位相 が 異な る も の (c1k gen45) 
Output [23:0] c1f ou: 
Wire pd1p, pd1n, fdp, fdn: 
Wire ckO0p, ckOn, ck45p, ck45n, ok90p, ck90n, ok135p, 
ck135n: 


wire [23:0] oc1f ou: 


// エ ュー エ 3 (pfd top, cp_1pE, vco chn) 用 パラ メー タ 
parameter TM PRCS=0.001: 


// [ns] , 時 間 精 度 


// エ 2 と r3 (cp_1pf, vco_chn) 用 パラ メー タ 
parameter K RTOT=524288.0: 
//rea1 変数 か ら nteger 変数 へ 変換 係数 


// エ 1 (pd top) 用 パラ メー タ 
parame も ter TD PD1=1.0: 
// [ns] , 遅延 時 間 , pp 
parame も ter TRE PD1=0.2: 
// [ns] , 立ち 上 が り ・ 立 ち 下 が り 時 間 , pp 用 
parameter K DITV PD1=50: 
// (TD+TREF) /K D+V 値 を 処理 サイ クル (cyc1e) と する . ED 用 
parameter TD FD=1.0:  //[ns], 遅延 時 間 , gp 用 
parameter TRE FD=0.2: 
// [ns] , 立ち 上 が り ・ 立 ち 下 が り 時 間 , FDp 用 
parame も ter K DTV EFD=50: 
// (TD+TRF) /K DTIV 値 を 処理 サイ クル (cyc1e) と する . FD 用 


// エ 2 (cp_1pf ) 用 パラ メー タ 
parameter CUR CP PD1=5.0: 
// [uA], チャ ー ジ ・ ポ ンプ の 電流 値 , pp 用 
parameter RTO B A PD1=1.0: 
// チ ャ ー ジ ・ ポ ンプ の 吸い 込み 電流 と 吐き 出し 電流 の 比率 
(Un-match の 影響 度 ) , pp 用 


parameter CUR CP ED=5.0: 
// [uA], チャ ー ジ ・ ポ ンプ の 電流 値 , rp 用 
parameter RTO _B A _FD=1.0: 

// チ ャ ー ジ ・ ポ ンプ の 吸い 込み 電流 と 吐き 出し 電流 の 比率 
(Un-match の 影響 度 ) , FD 用 


parameter OFS CUR CP _FD=0.0: 
// [uA] , 中 間 レ ベル 入力 で の チャ ー ジ ・』 電流 値 , FD 用 
parameter CP=8.0: 
// [pF], rpFr の キャ パシ タ 値 (フィルタ 定 数 ) 
parameter C2=0.2: 
// [pF], チャ ー ジ ・ ポ ンプ 切り 替り 時 の フィ ー ド スル ー 
緩和 用 キャ パシ タ 値 , rp で 使 


クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 


ディ ジタル ・ シ ミュ レー タ で 検証 する 


本 モジ ュー ル は デバ イス ・ ノ イズ を 考慮 すべ き ブ ロッ ク で 
す が , 今回 は 類 雑 さ を 避 ける 目的 か ら モ デリ ング を 省略 し 
て いま す . また , 図 3 a) 以 外 の 端子 と し て , リス ト 3 で 
は リセ ッ ト 信太 rst) と 電源 vdq, vss) 端子 を 追加 し て 
いま 最終 的 な レイ アウ ト 情報 に 配慮 し た ). 
3) pfd_top モ ジュ ー ル ( pfd_top.v) 

10 に 示す よう に, pfd_top モジ ュー ル は PD 回 路 と FD 
回 路 を 記述 し た pfd モジ ュー ル と , ディ ジタル 遅延 回 路 を 
記述 し た dgtl_dly モジ ュー ル か ら 構 成 さ れ ま す . た だ し , 


parameter REP=1000.0: 
// [Ohm] , rpF の 抵抗 値 ( フ ィ ル タ 定数 ) 
parameter TS _RK LPF=0.001: 
// [ns], ルン ゲ ク ッ タ 処理 の 周期 , rp 用 
parameter N STP TNT LPE=0.1: 
// [ns], ルン ゲ ク ッ タ 計算 の ステ ッ プ 幅 の 初期 値 
(cg 出力 の 更新 周期 の 初期 値 ) 


parameter RELERR LiPF=0.001: 
// 許 容 相対 エラ ー 項 (ルン ゲ ク ッ タ 計算 ) 
parameter ABSERR LiPF=0.00001: 


// 許 容 絶対 エラ ー 項 (ルン ゲ ク ッ タ 計算) 


//r3 (vco_chn) 用 パラ メー タ 
parameter CYCLE VCO=1 

//Vco で の 出力 更新 時 間 (=TM PRCS*CYCLE VCO) 
parameter THETA TNTT=0.0: 

// [rad], VCo で の 初期 位相 


defparam T1 .TM PRCS=TM PRCS: 
defparam T1.TD PD1=TD PD1: 
defparam T1.TRF PD1=TRE PD1: 
defparam 1.K DTV PD1=K DTV PD1 
defparam T1.TD FD=TD FD: 
defparam T1.TREF FD=TRE FD: 
defparam 1.K DTV FD=K DIV ED: 


defparam 12.TM PRCS=TM PRCS: 
defparam 12.K RTOT=K RTOT: 

defparam T2.CUR CP PD1=CUR CP PD1: 
defparam 12.RTO B A PD1=RTO B A PD1: 
defparam 1T2.CUR CP FD=CUR CP_FD: 
defparam 12.RTO B A FD=RTO B A FD: 
defparam 12.OFS CUR CP FD=OFS CUR CP_FD: 
defFparam 2 .CP=CP: 

defFparam T2 .C2=C2: 

defFparam T2 .RP=RP : 

defparam 12.TS RK=TS RK LPF: 
defparam 12.N STP TNT=N STP TNT LPF: 
defparam 12.RELERR=RELERR LPF: 
defparam 12.ABSERR=ABSERR TLPF: 


defparam 13 .TM PRCS=TM PRCS: 
defparam 13 .CYCLE VCO=CYCLE VCO: 
defparam 3 .THETA TNTT=THETA TNTT: 
defparam T3.K RTOT=K RTOT: 


//--- ネ ッ ト リ スト 
a881gn Ckgenp=Ck45D: 
a881gn ckgenn=ck45n: 
pfEd op 11 (rg, datp, datn, ckOp, ck0On, ck45p, ok45n, 
ck90p, ck90n, ck135p, ck135n, 
pd1p, pd1n, fdp, fdn, vVqad, Y88) : 
cp_1pf 12 (rst, pd1p, pd1n, fdp, fdn, oc1fE ou 上 , vYdd, 
Y88) : 
Yoco_chn T3 (rg, c1f ou 上 , ckOp, ckOn, ck45p, ck45n, 
ck90p, ck90n, ck135p, ck135n, Vdd, V88) : 
re 上 ime T4 (st, okk45p, ck45n, datp, datn, outp, outn, 
YQd, Y88) : 


endmodu1e 
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リス ト 3 vco_chn モジュール vco_chn.v) 


//VCoO チ ェ ー ン 回 路 a0=4 .6984://a0 で a4: 4 次 多項式 近似 の 係数 。 


//Revised on 01:23 05/09/14 a1=0 .2565: 
/ /Yin , You0, You45, Yout90, Vout135 :: signed- inteder, a2=-0.9054: 
4D1-OVe エ + 20b1 ヒ 上 -Cor@e a3=1 .6146: 
a4=-0.548: 
modu1e end 
Yoco_ohn (rs ,Y1n , OU 上 0D , Out 上 0n , Out 上 4 5p , Out 上 45n , Out 上 90p , Out 上 90n, 
ou 上 13 5p , Out 上 135n , VQd , YS8) : a1wayg8 后 (cyole) begin 
1nput て ら 8: FE (ras) begin 
nput [23:0] Yin: theta=THETA TNTT: 
npu VQd,V88: outOp = 1!b0: 
OutDpu モ 。 Ou ヒ OD , Ou 上 0n , Out 上 4 5D , OU 上 45n , Ou モ 90Dp, Out 上 90n, oOu モ 0n = “ou0p: 
ou 上 13 5p , ou 上 13 5n : out45p = 1'b0: 
red Cu ヒ OD, Out 上 0n, Ou 上 45D, Out 上 45n , Out 上 90p , Out 上 90n, Out45n = “Ou キ 45D: 
Ou 上 13 5Dp , ou 上 13 5n : out90p = 1'b0: 
Fed [31:0] tmp32Db 1,tmp32D 2: out90n = ^out90p: 
reg [23:0] Vin2, Vin24D 1mtH, Yin24b 1mtL: ou 上 135p = 1'b0: 
Out135n = “ou 上 135D: 
1ntegder CyC1G: end ease begin 
rea1 mp エ 1, mp r2, PT2, heta, d1 ヒ 上 : // シ ング ル エ ンド 入力 vin に 対す る リミッタ 
rea1 mp エ 3a, 上 mp エ 3b, mp 3c, mp r3d, tmDp エ 4a, 1E (Vin [23] ==1 !b0 ) 
tmp_r4b, tmp r4o, tmp r4d: Y1n2= (Vin > Yin24b 1mEH) ? Y1n24D 1mtH : vin: 
rea1 a0, a1, a2, a3, a4, Frq, Y1n 1: e1ge 


Y1n2=Y1n24D 1mtTi: 
parameter TM PRCS=0.001: 


// [ns], 時 間 単 位 // 整 数 か ら 実数 へ 変換 。 
parameter CYCLE VCO=1 : tmp32b 1={ {8{vin2[23] ) }), vin2) 
//Vco の 出力 更新 期間 を TM PRCs*CYCLE tmp r1=$1tor (tmp32b 1) : 
Vco で 与え る Yin エ 1=tmp r1/K RTOT: 


parameter THETA TNITT=0.0: 


// [rag] , 初期 位相 // 入 力 電圧 を , 周波 数 , 位相 へ 変換 
parameter K RTOT=524288.0: Frq = a0 
// 整 数 値 と 内 部 計算 用 実数 値 の 変換 比率 +a1*Y1n エ 1 
// 整 数 2^20 に 対し て , 実数 2.0 が 対応 する +a2* (Vim エ 1* マ 1n 1) 
nitia1 begin +a3* (Vnm 了 エ ]1* マ 1m エ ]1*V1n 1) 
PT2 = 6.28318530717958647692528676655901: +84* (V1mn 1*V1n 1*V1nm エ 1*V1n 1): 
theta=THETA TNTT: tmp_ エ 2=thea +dQ1 上 *PT2 ォ モエ ロ : 
d1 ヒ 上 =TM PRCS*CYCLE VCO: theta= (tmp_r2 > PT2)? (tmDp r2-PT2) : 上 mp エ 2: 


CyCle=CYCLE VCO: 
// 位 相 ( 四 つ の 位相 ) を 電圧 に 変換 し , ディ ジタル 化す る 


Yin24b 1mtH=24「!hOE 6666: tmp_r3a=$sin (theta) : 
Yin24D 1mtL=24'!h00 0000: tmp r3b=$gin (theta- (PT2/8.0) ) : 
tmp r3o=$gin (theta- (PT2/4.0) ) : 


リス ト 4 pfd_top モ ジュ ー ル pfd_top.v) 


// 位 相 ・ 周 波数 検出 回 路 + ディ ジタル 遅延 defparam 12.TM PRCS=TM PRCS: 
//Revised on 22:26 05/09/25 deEparam -TD=TD PD1 : 
//pd1 outp, pd1 outn, fd outp, FEd outn :: 1bit-data deFparam -TRF=TRE PD1 : 
defparam -K_DTV=K DIV PD1 
~1nc1ude "pFd.Y リ 
~1no1ude "dgt1 d1y.Y リ defparam -TM PRCS=TM PRCS: 
defparam -TD=TD PD1: 
modu1e pfd top (rst, datp, datn, okODp, ckOn, ck45p, defFparam -TRF=TRF PD1: 
ck45n, ck90p, ck90n, ck135p, ck135n, pd1 oup, defFparam -K_DTV=K DIV PD1 
pd1 outn, Ed ouEtp, fd outn, VQdd, YS8) : 
input rs, datp, datn, ckOp, ckOn, ck45p, ck45n, defparam T4 .TM PRCS=TM PRCS: 
ck90p, ck90n, ck135p, ck135n: defparam T4 .TD=TD FD: 
nput VQd, VS88: defparam T4 .TRF=TRE FD: 
Output pd1 Outp, pd1 outn, Fd_ouEp, Ed _outn: defEparam T4.K DIV=K DIV_FD: 
wire pd1 outp, Pd1 outn, pd1 ouEp1, pd1 ouEn1 
wire fd outp, fd outn, Ed ouEp1 , fd oun1 : defparam T5.TM PRCS=TM PRCS: 
defparam 1T5.TD=TD FD: 
//r2 て ェ 5 (dgt1 _d1y) 用 パラ メー タ defFparam T5.TRF=TRE FD: 
parameter TM PRCS=0.001: //[ngs], 時 間 単 位 deFparam T5.K DIV=K DIV FD: 


// エ 2 と +3 (two dgt1 q1y) 用 パラ メー タ // ネ ッ ト リ スト 

parameter TD PD1=0.2: // [ns] , 遅延 時 間 pfFd T1 (rg, datp, datn, ckOp, ckOn, ck45p, ck45n, 

parameter TREF PD1=0.2:  //[ns], 立ち 上 が り ・ 立 ち 下 が り 時 ck90p, ck90n, ck135p, ck135n, 

parameter K DTV PD1=50 : pd1 outp1 , pd1 ouEtn1, Fd ouEp1 , fd _outn1, Yad, VSS) : 
// (TD+TRF) /K DrV 値 を 処理 サイ クル (cyc1e) と する , PD dgt1 d1y 12 (pd1 oup1 , pd1 outD, vVdd , VS8) 

dgt1 d1y エ 3 (pd1 outn1 ,pd1 outn,Ydd, VS8) : 

// エ 4 と +5(two dgt1 q1y) 用 パラ メー タ dgt1 d1y T4 (Ed outp1 , Ed _outDp, Vdd, YS8) : 

parameter TD FD=0.2: // [ns] , 遅延 時 間 dgt1 d1y T5 (fd outn1,fFd outn,Ydd, YS8) : 

parameter TRF FD=0.2: // [ns] , 立ち 上 が り ・ 立 ち 下 が り 時 

parameter K DTV FD=50 : endmodu1e 
// (TD+TRF) /K DrV 値 を 処理 サイ クル (cyc1e) と する , ED 
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クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 


ディ ジタル ・ シ ミュ レー タ で 検証 する 


tmp_r3d=$g1n (theta- (PT2*3/8 .0) ) : 
tmp 4a=tmD 3a*K RTOT: 
tmp_r4D=tmD 3D*K RTOT : 
上 mp_ エ 4 で =tmD 3C*K RTOT: 
て tmp _ エ 4d=Emp エ 3Q*K RTOT: 
て tmp32b 2=$ エ to1 (tmp 4a) : 


ou て 0p = (tmp32b 2[31] ==1'b0)? 1!b1: 1!b0: 
Out0n = “ou0p: 

tmp32b 2=$ エ て o+1 (tmp 4b) : 

Out45p = (Emp32b 2[31] ==1'b0)? 1'b1: 1!b0: 
out45n = ~ て ou45D: 


tmp32b 2=$ エ て o1 (tmp r4c) : 
Out90p = (Emp32b 2[31] ==1'b0)? 1'b1: 1!'b0: 


out90n = ~out90D: Vc7f V] 較 
て tmp32b 2=$ ェ to1 (tmp r49) : 3 昌 志和 
out135p = (tmp32b 2[31] ==1'b0) ? 1'b1: 1'b0: 図 11 vco_chn.v モ ジュ ー ル に お ける 入力 電圧 と 出力 周波 数 の 関係 
out135n = ~out135p: Vco_chnv モ ジュ ー ル で は , 4 次 多項式 に よっ て 近似 し , 本 図 の 関係 を 実現 
Sn4 し て いる . 
end 
endmodu1e 


リス ト 5 pfd モジ ュー ル pfd.v) 


// 位 相 ・ 周 波数 検出 回 路 end 
//Revised on 03:11 05/09/25 end 
//pd1 outp, pd1 outn,Fd outp,Fd _outn :: 1bit-data endmodu1e 


modu1e // 位 相 検出 
pfd (rg , datp, datn , Cck0p, ck0n, ck45p, ck45n, ck90p, ck90n, modu1e pd(rg, da , Ok0 , Ck1,pd _ ou) : 
Ck135p , Ck135n,pd1 outp, pd1 outn, fd outDp,Fd outn, nput rg , da , ok0O,Ck1 : 


YQd , VS88) : 


ュ input て 5 , datp, datn, ckOp, ck0n, ck45p, ck45n, ck90p, ck90n, 


Cck135p, Ck135n : 
1nput YOd,V88: 
Outpu て pd1 outp, pd1 outn, fFd outp, Ed outn: 
wire pd1 outp, pd1 outn, FEd ouEtp, Fd outn: 
wire da , ck0, ok45, ok90 , ck135: 
wire Pd1 Oou: 
wire [1:0] fd ou1: 


ag881qn da= (datp==~datn) ? datDp : 1'b0: 

a881qn Ck0=(CkK0p==^ok0n)? ckOD : 1'b0: 

a881qgn Ck45= (ck45p==^ok45n) ? ck45p : 1'b0: 
a881qn Ck90=(ck90p==~ck90n) ? ok90p : 1'b0: 
a881qn CKk135=(Ck135p==^ok135n) ? ck135p : 1'D0: 
a8819n Pd1 OutDp=pd1 ou : // 差 動 正 側 
a881qn pd1 Outn=^pd1 outD: // 差 動 負 側 
aggtgn Fd outp=Fd ou 上 1 [1] : // 差 動 正 側 


agg1gn Fd outn=FfFd ou1[0] : // 差 動 負 側 


fd エ 1 (rg, da , Ck0 , ok45 , ck90 , ck135 , pd1 ou 上 ,FQ ou1) 』 
endmodu1e 


// 周 波数 検出 
modu1e fd(rgs, da , Ck0 , ok45 , ok90 , ck135 , pd1 ou,Fd ou) 』 
1nput 5 , dat,oCkO0,oCk45, ok90, ok135: 
output [1:0] fd ou: 
Output Pd1 ou』 
re [1:0] fd _ ou,tmp 2b: 
wire pd1 Out,pd2 ou: 


pd T1 (rg , da , ok45 , ck135,pd1 ou) : 
pd 12 (rg, da , ok0,ock90,pd2 _ ou) 


a1way8 @ (pd1 ou or posedge rst) begin 
せ f (エモ = ニ =1 ! D1 ) 
fd _ out <= 2'b00: 
e1se begin 
1 (pd1 ou ヒ ==1『「D1 ) 
て tmp_2D= (pd2 ou キモ ==1'b1) ? 2!b10 : 2!D11: 
e1ge 
tmp_2D= (pd2_ ou キモ ==1'b1) ? 2!b01 : 2!Db11: 
fd ou <= mp 2D: 


Output Dd _ ou: 
wire ok1,oCkq,pd ou 


detFF T1 (rs, da , Ck0 , ck1 ) : 

detFF 1T2 (rs, da , ok1 , Ckd) : 

detEfF ng 13 (rs, okq, ok1,pd ou) : 
endmodu1e 


// フ リッ プ フ ロ ッ プ 回 路 

modu1e deFf (rg , da , Ck , ou) : 
nput rg , da , CKk: 
Output の OU: 
ェ reg ou , ok1i , okdz 


a1way8 @ (ck or posedge rgt) begin 
せ F (エモ = ニ =1 !D1 ) 
QuE <= 1!b0: 
e1ge 
Ou 七 <= da 七 : 
end 
endmodu1e 


// 反 転 型 の フリ ッ プ 
フロ ッ プ 回 路 


modu1e detff ng (rg , da , Ck , Ou) : 
nput r ら 5 ,da ,Ck: 
Outpu の ou : 
Fed ou,Emp 1Db: 


a1way8 @ (ck or posedge rst) begin 
せ F (エモ = ニ =1 「D1 ) 
QuE <= 1!b0: 
e1se begin 
tmp_1D = (ok==1'b1)? da : ~da: 
Ou ヒ <= mp 1b: 
end 
end 
endmodu1e 
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リス ト 6 dgtl_dly モ ジュ ー ル ( dgtl_dly.v) 


// デ ィ ジ タル 遅延 回 路 
//Revised on 15:26 05/09/24, , , Digita1 De1ay 
//in, out :: 1b1E 


modu1e dgt1 d1y(in, ou ,Ydd, YS8) : 
npu VQd,V88: 


1nput 1n: 

Outpu の Ou: 

エ red Ou: 

reg Cn gr: //Counter 開始 フラ グ 
re [1 :50] 1n pr: //1 て 50c1ock 前 の in 値 


integer cycle, mp 11, tmp 12: 
nteger kk d1V2: 


nteger cn 上 drf, cntuD 上 エト : 
integer 1ndz: 
rea1 mp エ 1, mp エ 2: 


// [ns], 時 間 単 位 


dgtL_dy.v へ の 
入力 信号 図 


0.15X Vpp 


parameter TM PRCS=0.001: 


0.85X pp 
dgtlL dly.vy の b 
出力 信号 図 


| | 1ー ーー 

遅延 の い 遅延 い 
立ち 上 が り / 較 立ち 上 が り / 図 
立ち 下がり 凶 立ち 下がり 較 
時 間 7j。 時 間 7jF 


図 12 dgtl_dlyvy モ ジュ ー ル の 遅延 7)/ と 立ち 上 が り / 立 ち 下 が り 時 間 
7x』 の 定義 

dgt_dlyv モ ジュ ー ル の パラ メー タ の うち , 遅延 時 間 7。 と 立ち 上 が り / ち 
下がり 時 間 7 の 定義 を 行う . 今回 , 1 ビッ ト 期間 0.1ns よ り も 回 路 遅延 の 
ほう が 大 きい 場合 を 想定 し て , と 7 な の 区 別 を 行っ た . つま り , 7 な より 
短い パル ス は dgt_dly.vvy モ ジュ ー ル を 通過 で き な い . その た め , 図 中 の 記 と 
7 の 定義 は ユニ ー ク ( 一 意 的 ) に な っ て いる . 


リス ト 7 cp_lpf モ ジュ ー ル ( cp_lpf.v) 


// チ ャ ー ジ ・ ポ ンプ 回 路 キ ロー パス ・ フ ィ ル タ 回 路 
//Revised on 01:48 05/11/22 
//pd1p, pd1n, Edp, Edn, vou 


:: 符号 付き integer, 
4 ビッ ト ・ オ ー バ + 20 ビ ッ ト ・ コア 
modu1e cp 1]pf(rst, pd1p, pd1n, fdp, fdn, vou 上 , Ydd, 
Y88) : 

nput エ 8E: 
ュ input pd1p, Pd1n, 
nput YOd, Y88: 
Output [23:0] Vou 七 : 
reg rsE rk: // ル ン ゲ クッ タ 計 算 の 初期 化 フ ラグ 
ェ reg ca] rk: // ル ン ゲ クッ タ 計 算 の 有効 化 フ ラグ 
re OoOd_e エ エ =1!D0: 

// ル ン ゲ クッ タ 計 算 中 で エラ ー が 十分 小さ いか 否 か の フラ グ 
reqd Emp1b 1: 
red [23:0] Vou: 
re [31:0] hrs1 32D, 


fdp, fdn: 


thrs2 32b, tmp32D 1, 
て tmp32Db 1mtH, tmp32D 1mtT: 


ュ 1nted9er CyC1G: 
// ル ン ゲ クッ タ 処 理 の 実行 頻度 の 逆数 
integer h n, h n0: 


// ル ン ゲ クッ タ 計算 で の スキ ッ プ 数 , スキ ッ プ 数 の 初期 値 
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// [ns], 遅延 時 間 

// [ng] , 立ち 上 が り ・ 立 ち 下 が り 時 間 

// (TD+TREF) /K DrV 値 を 処理 サイ クル 
(cyc1e) と する 


parameter TD=0.2: 
parameter TREF=0.2: 
parame セ ter K DTV=50: 


ュ initia1 begin 
tmp_11=K _DTV: 
fk_div2=(tmp 11 > 50)? 50:tmp 1i1://L1mitter For K DIV 


tmp_ エ 1=(TD+TRF) / (k div2*TM PRCS) : 
tmDp_11=tmD エ 1: 
cyc1e=(tmp 1 < 1)? 1: tmp 11: 


tmp_ エ 1=CyC1e: 

tmp r2=TRF/ (tmp _r1*TM PRCS) : 
tmp_11=tmDp エ 2: 

cntup 上 rfF=(tmp 11 < 1)? 1: tmp 11: 


cn 上 QrF=0: 
Cnt g 上 エモ =1!D0: 
end 


a1way8 后 (cycle) begin 


For (ind=1 : ind<k d1V2 : 1nd=1nd+1 ) 
<= 1n pr[ind+1]:  //in の 遅延 
tn pr[k div2] <= in: 

Out <= 1n p エ [1]: 


in pr[indg] 


1F (cnt 8 ヒヒ ==1'b0) begin 
cn 上 上 QrfF=0: 
tmDp_11=(1n==1n pr[k d1v2] ) ? 0:1: 

end e1se begin 

Cn 上 上 QrF=cn 上 drf +1: 


FE (cn drf < cntup trf) begin 
EE (1 ロー=1nm pr[k d1vY2] ) 
tmDp 11=1 : 
e1se begin 
tmp_11=0: 
Eor (1nd=1 : ind<=onE 上 drf : nd=1nmd+ 1 ) 
//cntup_trf 未満 の 周期 を 持つ パル ス を 削除 
tn pr[k div2-1nd] <= 1n: 
end 
end el1ge 
tmp_11=0: 
end 
Cnt grE <= tmDp 11: 
end 


endmodu1e 


nteger cn dwn rk: 
// ル ン ゲ クッ タ 計算 に お ける スキ ッ プ 数 の カウ ンタ 
nteger hi min h: 
// ル ン ゲ クッ タ 計算 に お ける 最小 時 間 刻 み の 使用 頻度 
(cp 出力 更新 時 を 除く ) 


nteger Emp 11, mp 12: 

rea1 1Dp pd1, in pd1, ip fd, im fd, 

rea1 1 pd1, 1 fdz 

rea1] て] iin, a11, a12, a21, a22, b1: 

rea1 て] Db1 iin, て] x1, エ ] x2, k1a, k1b, k2a, k2Db, k3a, 
k3b: 


in fd: 


rea1 dif xx1, d1f x2, er x1, er エ xx2: 
エ @a] エ ] エ re]1er エ , エ ] abDSe エ エ : 
rea] mp エ 1, mDp 2, tmD エ 3, mp 4: 
parameter TM PRCS=0.001, // [ns], 時 間 単 位 
parameter K RTOT=524288.0: 
// 整 数 値 と 内 部 計算 用 実数 値 と 変換 比率 . 整数 2^20 に 対し て 
実数 2 .0 が 対応 
parameter CUR CP PD1=5 .0: 
// [nA] , 位相 検出 (pp) に 関わ る 充 放 電電 流 値 , 
具体 的 に は 正 側 電流 値 


parameter RTO B A PD1=1.0: 


BE 
// 負 側 出力 電流 値 の 正 側 出力 電流 値 に 対す る 比 (PD) 


parameter CUR CP EFD=5.0: 
// [uA] , 周波 数 検出 (FD) に 関わ る 充 放 電電 流 値 , 
具体 的 に は 正 側 電流 値 
parameter RTO _B A FD=1.0: 

/ / 負 側 出力 電流 値 の 正 側 出力 電流 値 に 対す る 比 (FD) 
parameter OFS CUR CP _FD=0 .0: 

// [uA] , 周波 数 検出 (FD) に 関わ る 中 間 レ ベル 電流 値 
parameter CP=8.0: // [pr], rpPF の キャ パシ タ 値 
parameter C2=0.2: 

// [pF], rpPF の キャ パシ タ 値 (クロ ッ ク ・ フ ィ ー ド スル ー 抑 制 用 ) 
parameter RP=1000.0: // [Ohm] , rp の 抵抗 値 
parameter TS _RK=0.001: // [ns], ルン ゲ ク ッ タ 処理 の 周期 
parameter N _STP 1TNT=0.1: 

// [ns], ルン ゲ ク ッ タ 計算 の ステ ッ プ 幅 の 初期 値 

(rpr 出力 の 更新 周期 ) 
parameter RELERR=0.001:  // 許 容 相対 エラ ー 項 
parameter ABSERR=0.00001: // 許 容 絶 対 エラ ー 項 


ュ initia1 begin 
エ ] re1err=RELERR : 
1] abSer エ =ABSERR : 
TS YK=1!D1 
hit minm h=0: 


// サ イク ル 数 の リセ ッ ト 
// 最 小 時 間 刻 み の 使 用 頻度 を 0 に セッ ト 


tmp エ 1=TS _RK/TM PRCS: 

tmp 11=tmDp エ 1: 

cyc1e=(tmp 11 < 1)? 1:tmp 11: 
tmp_r2=cCyc1e: 

a11=- (TM PRCS*1000 .0*tmp エ 2) / (C2*RP) : 


a12=-a11: 
a21= (TM PRCS*1000 .0*tmp 2) / (CP*RP) , 
a22=-a21: 


b1= (TM PRCS*1000.0*Emp r2) /C2: 


tmp_r1=cyc1e: 

tmp r2=(N_STP TNIT/(TM PRCS*oyoc1e) ) : 
tmp 11=tmp エ 2: 

h _n0=(tmp 11 < 1)? 1:tmp 11: 


1p_pd1=CUR _CP PD1*0.000001: 
// [A], PD に お ける 正 側 出力 電流 値 
1n pd1=-RTO_B A PD1*CUR CP PD1*0.000001: 
// [A] , PD に お ける 人 負 側 出力 電流 値 
1p_fFd=CUR CP_FD*0.000001: 
// [A], FD に お ける 正 側 出力 電流 値 
1n fd=-RTO_B A FD*CUR CP FD*0.000001: 
// [A] , FD に お ける 人 負 側 出力 電流 値 
1m fd=OFS CUR CP FD*0.000001: 
// [A] , FD に お ける 中 間 レ ベル 出力 電流 値 (理想 は ゼロ ) 
と] x1=0: 
エ ] x2=0: 


tmp32b 1mtH={8'h00, 24'h7F FEFFF) : 
tmp32b 1mt エ =~tmp32D 1mtH: 
end 


// チ ャ ー ジ ・ ポ ンプ 

a1way8 @ (pd1p or fdp) begin 
1 (pd1p==1!b1 ) 1 pd1=1Dp Pd1 : 
e1se 1 Dpd1=in Pd1 : 


3E((Fdp==1'b1) g(Edn==1'b0) ) 1 fd=ip fd: 
e1se 1F((Fdp==1'b0) g(Edn==1'b1 ) ) 1 Ed=in fdz: 
e1se 1 fd=im fd: 


// ル ン ゲ クッ タ 用 
1 11n=1 pd1-1 fFd: 
1 b1 11n ロ =b1* 了 エ ] 11n: 
TSt rk <= 1'Db1: 

end 


//Low Pasg F11te エ 
a1wayg 后 (cyole) begin 
FE (rS ヒ ==1'b1 ) begin 
1] x1=0: 
エ ] x2=0: 


// ル ン ゲ クッ タ 用 
ェ 1 b1 iin=0.0: 
rg TK <= 1!D1: 
h n=h n0: 

end e1se begin 


// ル ン ゲ クッ タ : 2 次 と 3 次 の 埋め 込み 式 


クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 


ディ ジタル ・ シ ミュ レー タ で 検証 する 


FE (rs 上 rk==1'b1) begin 


tmp 11=h n-onE dwn rk: 

tmp 12=(tmp 11 < 1)? 1: tmp 11: 
tmp_r1=h nz 

tmDp_ エ 2=tmp 12: 

tmp 3=tmDp 2/tmp 1: 

ェ 1 x1 = r] x1 +tmp r3x*k2a:  // 状 態 値 の 更新 
ェ 1 x2 = r1 x2 +tmp r3x*k2b:  // 状 態 値 の 更新 


h n=h n0: 
cnt dwn rk=h n0: 
Ca1 rk=1!D1 


end e1se begin 


FE (cn dwn rk==0) begin 


ェ ] x1 = エ ] xx1 +k2a: // 状 態 値 の 更新 
エ ] x2 = エ ] x2 +kK2D: // 状 態 値 の 更新 


hn = (hn << 2): 
Cnt dwn rk=h nz 
Ca1 rk=1!D1 : 

end e1ge 
Ca1 rk=1!b0: 


end 


FE (Ca]1 rk==1'b1) begin 


whi1e (oca] エ k==1'b1) begin 
// ル ン ゲ クッ タ 計算 の コア 部 
tmp_r3=h nz 
k1a= (a11* エ ] xx1 +a12* エ ] x2 + エ ] b1 11n)*EmD 3: 
k1D= (a21* エ ] x1 +a22* エ ] x2)* ヒ EmD 3: 


tmp_ エ 1 ユ = 了 1] x1 +k1a*0.5: 
tmp_ エ 2= 了 1] x2 +k1D*0.5: 
k2a= (a11*mD エ 1 +a12*EmD エ 2 
+ エ 1 b1 1in) *Emp エ 3: 
k2D= (a21*EmD エ 1 +a22*mD エ 2)*EmD エ 3: 


mDp_ エ 1 ユ 1= 了 1] x1 -kK1a +k2a*2 : 
mDp_ エ 2= 了 1] x2 -kK1D +k2D*2 : 
k3a= (a11*EmD エ 1 +a12*EmD エ 2 
+ エ 1 b1 1in) *Emp エ 3: 
k3D= (a21*EmD エ 1 +a22* ヒ mD エ 2)*EmD エ 3: 


dif x1=(-k1a +k2a*2.0 -k3a) /6.0: 
dif x2=(-k1b +k2b*2.0 -k3b) /6.0: 
tmp エ 3 = エ ]1 xx1 +(k1a +k2a*4 +k3a) /6.0: 
tmp r4 = 1 2 +(k1D +k2b*4 +k3D) /6 .0: 


tmp_ エ 1=(tmD エ 3 > 0)? tmD エ 3:-EmD エ 3: 
tmp_r2=(tmp r4 > 0)? tmDp エ 4:-EmDp 4: 
エエ 1= ヒ mD エ 1* 了 エ ] エ G1eG エ エ + 了 ] aDSer エ : 
ら エ エ xx2= ヒ mD エ 2 ォ エ ] エ G1e エ エ + 了 ] aDSer エ : 
tmp _ エ 1 ユ =(dif x1 > 0)? dif x1:-dif x1: 
tmp _ エ 2=(d1f x2 > 0)? dif x2:-dif x2: 


OOd _@e エ エニ = ( (エエ xx1 >= 上 mD て 1) g&& (err XX2 >= 
tmp_r2))? 1!b1:1!b0: 
if ( (good_err==1'b1) | | (h n==1) ) begin 
rst rk <= 1'b0: // 初 期 化 を 解除 
ca1 rk = 1'b0: // ル ン ゲ クッ タ 計 算 ル ー プ か ら 出る 
end e1se begin 
tmp 11 = hn >> 2: 
h n=(tmp 11 < 1)? 1: mp 11: 
hit min h=(tmp 11 < 1)? (hi min h+1 
hit min h: 


cn dwn rk=h n: 
end 
end 


end 
cnt dwn rk=ont dwn rk-1,// ダ ウン カウ ンタ の デ ク リ メン ト 


end 


て tmp_ エ 3= 了 1] x1*K RTOT: 
tmp32b 1=$rEo1 (tmDp 3) : 
3E (tmp32D 1[31] ==1'b0 ) 
Vout <= (tmp32Db 1 > tmp32b 1mtH) ? 


tmp32b 1mtH[23:0] tmp32b 1[23:0], // リ ミッ タ 


e1ge 
Vout <= (tmp32Db 1 < tmp32b 1mt) ? 


end 


モ tmp32b 1mtL[23:0] tmp32b 1 [23:0] , // リ ミッ タ 


endmodu1e 
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10 で は pfd_top モ ジュ ー ル pdf_topy) の 構成 要素 と し 
て 1 個 の dgtl_dly モジ ュー dgtl_dlyv) を 示し まし た が , 
実際 に は 図 4 に お ける 差 動 出力 で ある V ゅ と り ぁ に 対応 し た 
4 個 の dgtl_dly モジ ュー ル が 存在 し ます . リス ト 4 p1122) 
が ソー ス ・ コ ー ド で す . top モジ ュー ル か ら 受け 取っ た パ 
ラメ ー タ 値 は , dgtl_dly モジ ュー ル へ と 引き 渡さ れ ま す . 
4) pfd モ ジュール pfd.v) 

pfd モジ ュー ル は 図 4 を 実現 し て いま す . リス ト 5 p.123) 
が ソー ス ・ コ ー ド で す . 図 4 以 外 の 端子 と し て , リス ト 5 
で は , リセ ッ ト 信 坪 rst) と 電源 vdd, vss) を 追加 し て 
いま す . 
5) dgtl_dly モ ジュ ー ル ( dgtl_dly.v) 

dgtl_dly モジ ュー ル は , 図 4 の 入力 信 哨 入力 デー タ や 


生成 クロ ッ ク ) か ら 差 動 出力 で ある \j ヵ , Vp まで の 遅延 を 
実現 し て いま す . リス ト  p.124) が ソー ス ・ コ ー ド で す . 
四 つ の パラ メー タ 値 を pfd_top モジ ュー ル か ら 受け 取り ま 
す . TM PRCS は シミ ュ レ ーション 時 間 単 人 ~time 
gcale で 与え た の と 同じ 値 ), TD や TRF は 遅 延 定義 は 図 
12 p.124) の と お り . TRF は パル ス 幅 の 短い パル ス を 消滅 
させ る た め の パ ラメ ー タ ), K DTIV は 処理 時 間 の 間隔 で す 
( (TD+TRE) /K DIV が 処理 時 間 の 間隔 を 決める. なお, 実 
行 時間 の 短縮 と 遅延 処理 の 精度 は トレ ー ド オフ の 関係 に あ 
る . また , K DTIV の 上 限 値 を 50 と し た ). 
6) cp_lpf モ ジュ ー ル ( cp_lpfiv) 
図 6 お よび 図 7 を 実現 し ます . リス ト 7 が ソー ス ・ コ ー 
ド で す . 14 個 の パラ メー タ 値 を top モ ジュ ー ル か ら 受 け 取 


《 時 間 刻み 幅 可変 の 精度 2 と 精度 3 の ルン ゲ ク ッ タ 法 に よる 埋め 込み 型 アル ゴリ ズム 》 


常 微 分 方 程 式 の 初期 値 問題 初期 値 が 与え られ た 場合 の 常 微 分 方 程 
式 を 解く 問題 ) と し て , ルン ゲ ク ッ タ 法 が あり ます . ルン ゲ ク ッ タ 法 


は オイ ラー 肖 精度 1) か ら 派生 し た も の で , 精度 8 まで の アル ゴリ ズ 
ム が 紹介 され て いま ず ?②. な お , ここ で 精度 と は , 演算 結果 と 真 値 の 
エラ ー 最大 エラ ー) が 時 間 刻 み 幅 の 何 乗 に 比例 する か を 意味 し て い 
ます . それ ゆえ , 演算 エラ ー に つい て , 精度 1 で あれ ば 時 間 刻 み 幅 に 
比例 し , 精度 8 で あれ ば 時 間 刻 み 幅 の 8 乗 に 比例 し ます . 精度 が 高い 
ほう が 時 間 刻 み 幅 を 大 きく で きま す が , 1 回 当たり の 演算 量 が 増え 
ます . 
今回 , 計算 精度 の 維持 と シミ ュ レ ーション 時 間 の 短縮 と いう 目的 に 
お いて , 精度 2 と 精度 3 の ルン ゲ ク ッ タ 法 の 演算 結果 の 比較 に より , 
時 間 刻み 幅 を 可変 に する 解 流 アル ゴリ ズム) を 用 いま す . 詳細 な 根 

拠 の 解説 は 参考 文献 2) に 譲り ます が , ひと と お り 説明 し ます . 
まず , 本 文 の 図 7 を 状態 方 程 式 で 表し ます . 図 7 に 示す よう に 状態 
変数 を x+ と xs に すれ ば , 状態 方 程 式 は 式 B-1) に な り ま す . な お , 

が その まま 出力 電圧 に な り ま す . 
前 きま 
gZY| C2*Ap 。 C2*p 
1 

の * 9 


に は 
C2 


次 に , 時 間 刻 み 幅 を ヵ と すれ ば , 精度 2 の ルン ゲ ク ッ タ 法 の 推測 値 
え ょ 精度 3 の ルン ゲ ク ッ タ 法 の 推測 値 え .j」, お よび 精度 2 と 精度 3 の 
推測 値 の 差 は , 式 B-1) の 結果 を 受け て 式 B-2) よ うに な り ま す . 

人 0 ニ : ゴ / た ガオ 


人 デ の 


え 


1 


ー 和 を 
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_ gw POPCCCPPPPPPPPPPPPPPPPPPPEPPPEPPPPPPPP ( B-2) 
ー の 

_ の 

ー の XaYy+05e 

_ 

ー の ya 双 -wA1t2s が ek3 


ん 1 
ん 2 
ん た 3 
に , 上 式 B-2) について, 与え られ た 精度 を 満足 する よう に 時 
み 幅 を 変更 し な が ら , 演算 を 実行 し ます . すなわち 次 の 手順 


EP1< 3) に 従い ます . た だ し , 本 文 の 図 7 の 入力 , の 更新 時 は 
] 刻 み 幅 を 最小 幅 へ と 無 条 件 に セッ ト し ます . これ は , 推測 エラ ー 


が 大 きく な る こと が 予想 され る か ら で す . 


TEP1: 7 王 ヵ の と き の 台 を も と に , 7 ニム トカ の と き の え と 負 
を 求め る 

eSTEP2: 次 式 B-3) で 定義 する 目標 エラ ー g ヶ を 計算 する . た だ し , 
AE7EAA と 4gSEAA は , それ ぞ れ 許容 相対 エラ ー, 許容 絶対 エラ ー 
で ある . 


2 アニ 妥 * CLEAR だ 4 が SPA ーー…ーーーーーーーーーー… ( B-3) 


eSTEP3: 上 式 B-3) の 結果 e7 み と | ュー 名 を 比較 し て , 時 間 
刻み 幅 を 更新 する . つま り , 以下 の と お り で ある . 

1 ュー 名 ミ 7 の と き , 77+1 の と き の 双 値 と し て 
STEP1 で 求め た 史 .、」 で 良し と する . 時 間 刻 み 幅 を ヵ つ 4X ヵ に 変更 
し た うえ で STEP1 を 実行 する 2 を 求め る ) 

1 ュー 宮 | っ g の と き , 時 間 刻 み 幅 を ヵ っ 025X ヵ に 変更 
し た うえ で STEP1 に 戻っ て 再度 え .」 を 求め る . た だ し , 時 間 刻 み 
呈 の 変更 値 が 最小 刻み 幅 と な っ た 場合 は , STEP1 で 求め た 兄 .、」 で 
良し と する . その うえ で , 時 間 刻 み 幅 は その まま で STEP1 を 実行 
する (jr を 求め る ). 


り ま す . TM PRCS は シミ ュ レ ーション 時 間 単 侯 ~time 
scale で 与え た の と 同じ 値 ), K_ RTOT は 内 部 変数 実数 ) 
と 信 豆 整数 表現 ) の 比 係数 ) で す . CUR_CP PD1, 
CUR _CP EFD, RTO B A PD1, RTO B A FD, 
OFS CUR CP ED は , PD 用 も し く は FD 用 で , 素子 の ミ 
スマ ッ チ や オフ セッ ト に よっ で 図 6), 正 側 電流 や 負 側 電 
流 , も し く は その 中 間 レ ベル 電流 の 間 で 比例 関係 が 崩れ た 
場合 を 想定 し た パラ メー タ で す . それ ぞ れ の 出力 電流 の 算 
出 式 は , リス ト 7 の initia1 ブロ ッ ク 中 の ip pd1, 
in pd1, ip fd, in fd, im fd の 算出 箇所 に 相当 し ま 
す . また , CP, C2, RE は 図 7 の と お り で す . TS RK, 
N_8TP TNIT, RELERR, ABSERR は 図 7 の フィ ル タ 出 力 
( 常 微 分 方 程 式 ) を 数 値 計算 で 解く < 際 の パラ メー タ で ず 左 
掲 の コラ ム 時 間 刻 み 幅 可変 の 精度 2 と 精度 3 の ルン ゲ ク ッ 
タ 法 に よる 埋め 込み 型 ア ル ゴ リ ズム 」 を 参照 なお, 図 6 
お よび 図 7 以 外 の 端子 と し し て , リス ト 7 で は , リセ ッ ト 信 
号 rst) と 電源 vdd, vss) を 追加 し て いま す . 
7) re_time モジ ュー ル re_time.v) 

re_ time モ ジュ ー ル は , 図 10 の よう に vco_chn モ ジュ ー 
ル で 生成 し た クロ ッ ク に よっ て 入力 デー タ を 取り 込み ます . 
具体 的 な 回 路 は , 図 4 a) の PD 回 路 に 使用 し て いる 差 動 型 
フリ ッ プ フロ ッ プ 回 路 と 同じ こ に なり ます . また , リス ト 8 
が ソー ス ・ コ ー ド で す . 


リス ト 8 re time モジ ュー ル ( re_time.v) 


// デ ー タ 再生 成 回 路 

//Revised on 18:40 05/08/11 

//rst , ckp, ckn,dat inp,dat nn,dat outp,dat outn :: 
1b1 上 


modu1e re ime (rg , CkD , Ckn , dat 1nDp,da inn,dat ou て p, 
dat outn,YQdd,VS8) 』 
1nput rs , okp, okn , dat 1np,dat 1nn: 
1nput YOd,V88: 
Outpu て da oup,Qdat outn: 
reg da ou』 
wire da outD, dat outn, ck,dat in: 


a881qn da ヒ ou て tp= (dat ou モ ==1'b0)? 1'b0 : 1'D1: 
a881qn da Outn=~^da OutD: 


a881gn ck= (ckp==~ckn) ? ckp : 1!b0: 
a881qn da 1n= (dat 1np==~dat 1nn)? 
dat inp : da inn: 


a1way8 @ (posedge ck or negedge ck) begin 
ff (8 上 ==1 ! D1 ) 
dat ou <= 1!b0: 
e1ge 
dat out <= dat in: 
end 


endmodu1e 


クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 


ディ ジタル ・ シ ミュ レー タ で 検証 する 


⑯ 入力 デー タ 生 成 回 路 : 23 ビッ ト の 疑似 乱数 を 発生 

図 10 に お ける gat_genv は , 23 ビ ッ ト ( 229- 1 の 乱数 ) 
の 疑似 乱数 生成 回 路 で す . 疑似 乱数 は M 系 列 に よっ て 生成 
し て いま p.119 の コラ M 系 列 に よる 疑似 乱数 の 生成 」 
を 参照) リスト 9 が ソー ス ・ コ ー ド で す . 


人 測定 回 路 : リフ ァ レ ンス 信号 と 位相 を 比較 

10 に お ける msjty は 測定 回 路 で す . リス ト 10 が ソー 
ス ・ コ ー ド で す . strt 信号 が 'H' で 測定 状態 と な り , 入 
カカ デー タ 列 を 作る リフ ァ レ ンス 信号 reF clk と voo chn 
モジ ュー ル で 内 部 生成 し た クロ ッ ク gen c1k m90 の 位相 
を 比較 する こと で (理想 的 に は , reE c1k と 
gen clk m90 は 90 位相 が 異な る ), ジッ タ 量 を 測定 し ま 
ず 図 9). また , 5 個 の パラ メー タ 値 を top モジ ュー ル か ら 
受け 取り ます . TM PRCS は シミ ュ レ ーション 時 間 単 位 
( ~timescale で 与え た の と 同じ 値 ), CYCLE MT は 測定 
を 行う 際 の スキ ッ プ 数 , T AVRG は ジッ タ 量 の 平均 値 を 求 
め る た め の サ ンプ ル 期 間 , T AVMT は ピー ク ・ ジ ッ タ 量 や 
RMS ジ ッ タ 量 を 求め る た め の サ ンプ ル 期 間 , TS CLK は 入 
カカ デー タ を 生成 する た め の リ ファ レン ス ・ ク ロッ ク 周 期 で 
す . 


リス ト 9 疑似 乱数 生成 回 路 dat_gen.v) 


// 入 力 デー タ 生成 回 路 

//Revised on 11:37 05/08/20, , Data Gene エ aCr, 
//23b1it -M- ser1eg 

//pEd_outp,pFd outn :: signed- inmteger , 

// 4b1t-over + 20b1 上 -cCore 


modu1e da gen(rst,C1k, datp , datn) : 
nput rs ,C1k: 
Outpu セ datp , datn : 
reg XX PrDbS [1:23] : 
wire datp, datn , da : 


a881qn da ヒ =x prbs [23] : 
a881qn daD=da : 
a881qn da ヒロ =~da 七 : 


nteger 1ndex: 


a1way8 @ (oc1k or posedge rst) begin 
FE (rS ヒ ==1'b1 ) begin 
For ( index=1 : 1ndex<24 : 1ndex=1ndex+1 ) 
x prbg [index] <= 1'D1: 
end e1se begin 
For (index=1 : hndex<23 : 1ndex=1ndex+1) beg1n 
x prbs [24-1ndex] <= x prbs[23-index] : 
end 
メ prbg [1] <= (x prbg [18] 
end 
end 


ス 


メ DrDbg [23] ) : 


endmodu1e 
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Mo 
> 


SS 


リス ト 10 測定 回 路 msjt.v) 


// 測 定 回 


路 


//Revised on 14:19 05/09/03 

//Eva1uate peak to peak, rms ]1EEer between reF-o]ook 
and 135deg(45deg+90deg) -Yco-geno1k 

//But we assume て tha peak ]1EEer 1g wihn +/- 50psec(- 
90deg~+90deg) 


modu1e ms]t(strt,ref oc1k,gen c1k m90): 
nput strt,refF co1k,qen o1k m90: //gen cl1k m90 => Vou45 
Fed ref c1K pr,den Cl1k m90p エ : 
red en m] sh: 


if(cnt av < sg av) begin // 発 人 
tmp 11=F1g set gen oc1k m90: 
1E (tmp 11==0) begin 

m] cn=0: 
FE1g_ set gen cl1k m90=(den c 


k_m90 


本 den c1k m90p エ )? 0:1: 


end e1se begin 
m] Cnt=m] Cn+1 : 
FE1g set gen Cl1k m90=(ref c 


E((ref oc1k != refF o1k rr) 


k==refF C1K て )? 
R0 有 H 
gg& (m] cnt <= 
ts_180deg) ) begin 


1nteger 
1ntegder 
nteger 
1nteger 


CyC16: 

ts_av,Cn avY,g avm],tm avm]j,ts 90deg,ts 180deg: 
m] er エ : 

FE1g_set gen co1k m90: 


//Vco 生成 クロ ッ ク の エッ ジ で 1 と な り , リフ ァ レ ンス ・ ク ロッ ク 


integer 


の エッ ジ で 0 と な る 
FE1g _set gen c1k m90Esh: 


//VcCo 生成 クロ ッ ク の エッ ジ で 1 と な り , それ 以外 の Vco 状態 で 


1nteger 
ntegder 
ntegder 
1nteger 
ntegder 
1nteger 
ntegder 


0 と な る 
Cnt ry av 
m]_Cnt ,m] av エ , エ S1 av: 

tsh min gen c1k,tsh max gen c1k: 
cnt gh 
m] pk min,m] pk max: 
tmp_11,tmp 12: 

hit rms,m] rms2 1: 


ェ ea1 mp エ 1,mDp エ 2,EmD エ 3: 


parameter TM PRCS=0.001, //[ns], 時 間 単 位 
parameter CYCLE MJ=1: 

/ / 測 定 を 行う 際 の スキ ッ プ 数 . 1/CycLg Mr が 測定 頻度 を 意味 する 
parameter T AVRG=500: 

// [ns], strt が ->H 後 の ジッ タ 平 均 値 を 求め る た め の 時 間 を 指す 
parameter TT AVMJ=3000: 

// [Ins], ピーク ジッタ や RMs ジ ッ タ 値 を 求め る た め の 時 間 を 指す 
parameter TS CLK=0.2: //[ns], リフ ァ レ ンス ・ ク ロッ ク の 周期 


nitia1 begin 
gen c1k m90pr=1'b0: //1 c1ockk 前 の VCO C1ock (vout45) 
refF C1k pr=1'D0: //1 c1ock 前 の ReE C1ock 
en m] tsh=1'b0 : //Vco の 半周 期 測定 フラ グ を ogg に 設定 
(strt="H" 時 に ON と する ) 
CyC1e=CYCLE MO: 


tmp_ エ 1=T AVRG/ (TM PRCS*oyo1e) : 
ts_av=tmp r1: // 発 生 ジ ッ タ 量 の 平均 値 を 求め る 期間 (int) 
tmp_ エ 1=T AVMJ/ (TM PRCS*oyo1e) : 
ts_avm]j=tmp 1, // ピ ビー ク ・ ジ ッ タ や RMs ジ ッ タ 量 を 

測定 する 期間 (int) 
tm _avm]=ts_avm]j+ts av: //rms & peak ]1tter の 


算出 終了 時 刻 (int) 


F1g set gen Cc1k m90=0: 
F1g set gen c1k m90tsh=0 : 


tmp r1=TS _CLK*0.25/ (TM PRCS*oyoc1e) : 
tg_90deg=tmp_ 了 1 //90deg に 相当 する 期間 (int) 
ts_180deg=ts_90deg*2, //180deg に 相当 する 期間 (1nt) 


cnt try_av=0: //VCO 生 成 ク ロッ ク と リフ ァ レ ンス ・ ク ロッ ク の 
エッ ジ が 180geg 期 間 で ミー ト し た 回 数 を カウ ント 
tsh min gen clk=ts 180deg: //VCo 周 期 の 最小 値 (int) 


tsh max gen clk=ts 180deg: //VCOo 周 期 の 最大 値 (int) 


m] aV エ =m] avr +(m] cnt-s 90deg) : 
Cm ヒエ aV=Cn ヒ 上 エ Y avV+1: 
end 
end 
Cm 上 aV=CnE av+1 : 


end e1se 1F(cnE av == sg ay) begin 
// 発 生 ジ ッ タ 量 の 平均 値 を 計算 
1f (Cm 上 上 ry avV==0) begin 
Sdisp1ay ("Did not 1ook... 8o quit now.") : 
ド に 』 ス <) > 
end e1se begin 
tmD 1=m] avV エ : 
上 mD_ エ 2=Cn 上 上 エ Y avVj 
tmp エ 3=tmp エ 1/tmDp 2: 
8S1 七 aV=mD エ 3: 
end 
Cm 上 aV=CnE av+1 


end e1se 1iF(onE av > tm avm]) // 測 定 し な い 


Cnt aV=tm avm]+2: 


e1se begin  //"onE avV > Ss _ av" && "Cmt aV <= 
tm_avm] "を 満た す ケ ー ス 発生 ジッ タ の ピー ク 値 や RMs 値 を 測定 
//Vco で 生成 し た クロ ッ ク の 最小 半周 期 と 最大 
半周 期 を 測定 
tmp 11=F1g set gen c1k m90: 
FE (tmp 11==0) begin 
m] cn=0: 
FE1g set gen c1k m90=(gen cl1k m90 
==9en C1k m90pr) ? 0:1: 
end e1se begin 
m] Cnt=m] Cn+1 : 
FE1g_ set gen Cl1K m90=(ref C1k==refF C1K Fr)? 
0』: 
1F((ref c1k != ref oc1fkk pr) && (m] cnt <= 
ts_180deg) ) begin 
m]〕 e エ エ =m] Cm -s 90deg -rg1t av: 
m] rm82 1=m] rm82 1 +(m]〕 er エエ *m] @ エ エ ): 
m] pk m1n=(m] er エ < m] pk min)? m]〕 er エ 
m]_ pk min: 
m]〕 pk max= (m] er エ > m] pk max)? m]〕 er エ 
m] pk max: 


h1t 上 エ rms=h1 ] 上 エ ms+1: 
end 
end 


FE1g set gen c1k m90Esh=(gen c1k m90 
==9en Cl1K m90pr) ? 0:1: 
FE(F1g _ set gen C1kK m90 上 sh==0 ) 
cnt tsh=(en m] tsh==1'b1)? (cnt tsh+1) : 1: 
//Vco で 生成 し た クロ ッ ク 周 期 を カウ ント 


e1se begin 


Cnt aV=0: 
cn sh=0: 


m] Cn=0: 
m] av エ =0 : 


// 観 測 用 の カウ ンタ (#cycte ご と に 更新 ) 
//VcCO-C1ock エッ ジ か ら の 経過 時 間 を 対象 と し た 
カウ ンタ 
//jitter を 求め る た め の カ ウン タ 
//jitter の 平均 を 求め る た め の 総 和 (int) 


3F(en m]〕 tsh==1'b1) begin 
//Vco で 生成 し た クロ ッ ク 


tsh mn gen oc1k=(ocnE tsh < 


cn gh : 


周期 の 最小 と 最大 を 更新 
tsh min gen co]1k) ? 
tsh min gen o1k: 


m]〕 rms2 1=0, //]1tter の rms^2 を 求め る た め の 変 数 (int) 
m〕 pk _min=0,: //]1tter の peak 最 小 値 を 求め る た め の 変 数 (int) 
m〕 pk max=0, //]1tter の peak 最 大 値 を 求め る た め の 変 数 (int) 
hit trms=0, //]itter の rms^2 計算 に 使用 し た サン プル 数 
( in) 
end 


a1way8 后 (cycle) begin 
せ FF (8 ヒエ ヒ = ニ =1 !b0 ) cn av=0: 
e1se begin 
ref C1K pr <= ref Cc1k: 
gen cC1k m90pr <= gen cl1k m90: 
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tsh max den C1k=(cn tsh > sh max den c1k) ? 
cnt tsh : sh maxx den Cl1k: 
end e1gse 
en m] tsh=1!D1 : 
Cnt 上 sh=1 
end 
Cm 上 aV=CnE av+1 : 
end 
end 
end 


endmodu1e 


』 クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 
胃 。。) ティ ジタル ・ シ ミュ レー タ で 検証 する 


も 
志 
志 


図 13 

リス ト 1 に よる シミ ュ レ ーション 結果 

( 0 て 1.9s) 

vin_r( VCO 入 力 電圧 ) と fr VCO 発振 周 波 
数 ) が , 1.0s 付近 で ロッ ク し て る こと が わ 


か る . 


坊 


. 1 
。 46061219 


いり 


tb 64h4013E 


図 14 


リス ト 1 に よる シミ ュ レ ーション 結果 
( 1.0 ws 付近 ) 


datp と ck135p が 同期 し て いる の を 確認 で きる . 


『 シミ ュ レ ーション 結果 


今回 デバ イス ・ ノ イズ は 考慮 し て いま せん が , ディ ジ 
タル 遅延 値 や ロー パス ・ フ ィ ル タ の 定数 , チャ ー ジ ・ ポ ン 
プ の 電流 値 の 調整 し だ いで , ジッ タ 量 や ロッ ク す る まで の 
時 間 が 変わ り ま す . ここ で は , ディ ジタル 遅延 値 に 着目 し 
て , ジッ タ 量 や ロッ ク す る まで の 時 間 を 確認 し ます . な お , 
シミ ュ レ ーション 時 間 に つ いて は , 24GHz 動作 の Celeron, 
736M バ イト の メイ ン ・ メ モリ , Windows XR Home 
Edition,。 SP2) で リス ト 1 の テス ト ベン チ を VeritakWin で 
実行 し た と ころ , 7~ 8 分 で 処理 が 終了 し まし た . 


⑯ シミ ュ レ ーション し て 結果 を 確認 

リス ト 1 を VeritakWin で シミ ュ レ ーション し まし た . 
rs リセ ッ ト 信号 ), datW 入力 デー タ ), ck135% VCO 
で 内 部 生成 し た クロ ッ ク ), vin rT VOO 入 力 電 由 V], 
実数 値 ) Er VCO 発振 周 渡 数 GHz], 実数 値 )) お よび 
gtrt ms] 測定 開始 フラ グ ) に 着目 し て プロ ッ ト し た の 


VCO 入 力 電圧 [ V] 較 
1.0 6 


0.8 


0.6 


0.4 


0.2 


0 0.5 1.0 1.5 2.0 
時 間 [ ms]【 


図 15 リス ト 1 の パラ メー タ 値 で の シミ ュ レ ーション 結果 

リス ト 【 テス ト ベン チ ) の パラ メー タ 値 を 用 いた シミ ュ レ ーション 結果 . VCO 

入力 電圧 が ロッ ク す る まで の 時 間 ロッ ク 時 間 ) や ロッ ク 状態 に お ける リプ ル を 
認 で きる . 


が 図 13 と 図 14 で す . 13 は 0 て 19xs 期 間 の プロ ッ ト で 
す . vin r1 と frq は アナ ログ 波形 表示 で , 10s 付近 で 
すでに ロッ ク し て いる こと が わか り ま す . また , 図 14 は 10 
ん s 時 あたり を 拡大 し た も の で , gatp と ck135p が 同期 し 
て いる こと を 確認 で きま す . 

15 は , 外部 ファ イル rslt.dat へ 書き 込ん だ 結果 で す . 
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VCO 入 力 電圧 [ V] 図 
1.0 と 
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図 16 ディ ジタル 遅延 値 TD PD1 と TD_FD) を 2 倍 に し た 場合 の シ 
ミュ レー ショ ン 結 果 


VCO 入 力 電 上 中 図 13 の vin r1 に 相当 ) の 時 間 応 答 を 
gnuplot で プロ ッ ト し まし た . さら に , rslt.dat フ ァイル の 
最後 尾 に 測定 し た ジッ タ 量 が 書き 込ま れ ま ず 表 1 を 参照 . 
た だ し , 単位 は RMS の 2 乗 値 が ps2 で , これ と 許容 誤差 を 
満足 し な い 回 数 を 除い た 結果 の 単位 は ps). ジッ タ 量 の ピ 
ー ク 値 と RMS 値 を 表 2 に 示し ます . 目標 と し た ジッ タ 量 
を クリ ア し て いま す . 


⑯ ディ ジタル 遅延 値 を 変更 し て シミ ュ レ ーション 

リス ト 1 で 指定 し た パラ メー タ 値 の うち , ディ ジタル 遅延 
に 関す る TD_PD1 と TD ED の 二 つ の パラ メー タ を , それ ぞ 
れ 06 と 10 に 変更 し , シミ ュ レ ーション を 実行 し まし た . 

上 記 と 同じ よう に , VCO 入 力 電圧 の 応答 を 図 16 に 示し 
ます . 図 15 と 比べ て , ロッ ク し た 状態 の リ プル が 大 きい 
こと を 確認 で きま す . また , ジッ タ 量 の ビ ピーク 値 と RMS 
値 を 表 2 に 示し ます . 上 記 と 比べ て ジッ タ 量 は 大 きく なり , 
RMS 値 で 目標 を クリ ア で き て いま せん . 
今回 の 例 で は 燃 雑 さ を 避け る た め , デバ イス ・ ノ イズ を 
無視 し まし た . これ を 考慮 する 場合 に は , あら か じ め ノ イ 
ズ 信号 を 外部 ファ イル と し て 用 意 し て , 感度 の 高い VCO 
入力 に Sreadmemnm コマ ンド な ど で 加 える 方 法 が あり ます . 

また , シミ ュ レ ーション 時 間 を 短縮 する た め に , 時 間 単 
位 に 1ps を 指定 し まし た . 今 , リス ト 1 の テス ト ベン チ に 
て , 時 間 単 位 を 小さく , 例え ば 100fs に する と , RMS 値 
( ジッ タ 量 ) は 05ps 程度 に 下がり , ロッ ク に 入る 付近 の 大 
き な 振 動 が ほぼ ば なくなり まず ポ 図 15 を 参照 . 厳密 に は VCO 
の デバ イス ・ ノ イズ な どの 影響 も 考慮 する 必要 が ある ). ち 
な み に , 時 間 単 位 を 100fs と し た 場合 の テス ト ベン ヂ ( モ 


【 


All 


130 Design Wave Magazine 2006 November 


2.0 
時 間 [ ms]| 


表 1 測定 し た ジッ タ 量 


正 側 ピ ー ク 1 2 
負 側 ピー ク 1 ー 1 
RMS の 2 乗 値 0578716 
最小 値 99 
最大 値 101 


VCO 生成 クロ ッ ク 半周 期 


ルン ゲ ク ッ タ 計算 
で 許容 誤差 を 満足 0 
し な か っ た 回 数 


表 2 シミ ュ レ ーション 結果 ジッ タ 量 ) 
ディ ジタル 遅延 値 
( TD_PD1 と TD_FD) 備 考 
TD_PD1= 03ns リス ト 1 の パラ メー タ 値 を 採 
TD_FD= 05ns 用 . 図 15 も 参照 
左記 の ディ ジタル 遅延 値 以 外 
は , リス ト 1 の パラ メー タ 値 


TD_PD1= 06ns 
TD_FD= 1.0ns 9 を 採用 . 
図 16 も 参照 


ジュ ー ル 名 は tb_top100fs) を リス ト 11 に 示し ます . シ 
ミュ レー ショ ン 時 間 は 7 8 倍 に な り ま す . 興味 が あれ ば , 
シミ ュ レ ーション し て みて いた だ きた いと 思い ます . 
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あり は ら ・ え いい ち 
アナ ログ / ミ ックス ト ・ シグナル IC コ ン サ ル タ ン ト 


筆者 プロ フィ ー ル > 
在原 栄一 . アナ ログ / ミ ックス ト ・ シグナル |IC 設 計 や 機能 検証 に 関す 
る 受託 設計 や 技術 コン サル ティ ング に 従事 し て いる . 


クロ ッ ク ・ デ ー タ ・ リ カバ リ 回 路 の 機能 を 


ディ ジタル ・ シ ミュ レー タ で 検証 する 


リス ト 11 時 間 単 位 を 100fs と し た 場合 の テス ト ベン G( tb_top100fs.v) 


//CDR 回 路 用 テス ト ベン チ defparam T2.RTO_B A PD1=RTO B A PD1: 

//Revised on 01:29 05/11/23 defparam T2.CUR CP FD=CUR CP_FD: 
defparam 12.RTO_B A FD=RTO B A _FD: 

~timesca1e 100fs / 100fs //<== 時 間 単 位 100fs の 場合 に 変更 ! defparam T2.OFS CUR CP FD=OFS CUR CP_FD: 

defFparam T2 .CP=CP: 

~1nc1ude "EoOD .Y リ defFparam T2 .C2=C2: 

~1no1ude "dat gen.Y リ defFparam T2 .RP=RP: 

~1no1ude "ms] 上 .Y リ defparam T2 .CYCLE VCO=CYCLE VCO: 

defparam 12.THETA TNTT=THETA TNTT: 


modu1e 上 D Eop100Fg: 
reg TS 上 , ref clk: defparam 13 .TM PRCS=TM PRCS: 
red 8SrE ms] 上 , okgen: defparam 13 .CYCLE MJ=CYCLE MJ: 
red ms] oc1k: defparam T3 .T AVRG=T AVRG: 
re ms] 上 oc1k2: defparam T3 .T AVMJ=T AVMJ: 
wire Out D, ou n, okgen45 p, ckgen45 n, dat p, da nz defparam 13.TS CLK=TS _CLK: 
Wire VQQ, Y88: 
Wire [23:0] c1F ou //--- ネ ッ ト リ スト 

a881gn vVdd=1!D1 : 

nteger 上 S refFolk: a881gn VS8=1!Db0: 

integer mod: a881qn Ckgen45= (ckgen45 p==ckgen45 n)? 1'b0 

nteger mp 11,EmD 12: ckgen45 D: 

rea1 mp 了 1,EmD エ 2: dat_ gen 11 (rst,ref cl1k,dat p,dat n): 

top 12 (rst, da p, da n, ou p, ou n, ckgen45 p, 

// ト ッ プ 回 路 (top) と 測定 回 路 (msjt ) 用 パラ メー タ ckgen45 n, c1fF out, vqd, vV88) , 

parameter TM PRCS=0.0001: mg 上 3 (str ms]t, ref o1k, okkgen45) : 

// [ng] , 時 間 単 位 <== 時 間 単 位 100Es の 場合 に 変更 ! 


ュ initia1 begin 
// ト ッ プ 回 路 (top) 用 パラ メー タ tmp エ 1=TS CLK/TM PRCS: 
parameter K RTOT=524288.0: //rea1 変数 か ら integer 変 数 へ 変換 係数 ts _refo1k=tmp エ 1: 
parameter TD PD1=0.3: //[ns] , 遅延 時 間 , pp 用 // 入 力 デー タ 生成 用 リフ ァ レ ンス ・ ク ロッ ク 有 周期 
parameter TRF PD1=0.1: //[ns], 立ち 上 が り ・ 立 ち 下 が り 時 間 , pp 用 mcd=$Fopen ("rs1t .dat") : 
parameter K DTV PD1=20: // 結 果 書 き 込み 用 ファ イル の 指定 
// (TD+TREF) /K DrV 値 を 処理 サイ クル (cyc1e) と する . ED 用 て 8 キ =1「D1 // リ セッ ト 有効 
parameter TD EFD=0.5: //[ns], 遅延 時 間 , rp 用 8t エ t mg]=1'b0: // 測 定 無 効 
parameter TRF FD=0.1:  //[ns], 立ち 上 が り ・ 立 ち 下 が り 時 間 , rp 用 
parameter K DTV FD=20 : #1000 rs=1!D0 : 
// (TD+TRF) /K DrV 値 を 処理 サイ クル (cyc1e) と する . gp 用 // リ セッ ト 解除 <== 時 間 単 位 100Es の 場合 に 変更 ! 
parameter CUR CP PD1=125.0, //[uA], チャ ー ジ ・ ポ ンプ の 電流 値 , #12000000 strt mg]=1'D1 : 
PD1 // 測 定 有効 <== 時 間 単 位 1r00Fs の 場合 に 変更 ! 
parameter RTO_B A PD1=1.0, // チ ャ ー ジ ・ ポ ンプ の 吸い 込み 電流 #8100000 begin // <== 時 間 単 位 100s の 場合 に 変更 ! 
と 吐き 出し 電流 の 比率 (Un-match の 影響 度 ), pp 用 $Edisp1ay (mcd, "ジッ タ 量 -- - 平 均 値 : sq", エ 3 .rs1t ay) : 
parameter CUR CP FD=125.0, //[uA], チャ ー ジ ・ ポ ンプ の 電流 値 , gp 用 $Fdisp1ay (mcd, "ジッ タ 量 -- - 正 側 ピ ー ク 値 : sq, 
parameter RTO B A FD=1.0, // チ ャ ー ジ ・ ポ ンプ の 吸い 込み 電流 負 側 ピー ク 値 : *d", エ 3 .m]_pk max, 13.m〕 pk min) : 
と 吐き 出し 電流 の 比率 (Un-match の 影響 度 ), rp 用 tmp_ エ 1=T13 .m] rms2 1: 
parameter OFS CUR CP FD=0.0: tmp_ エ 2=T3 .hit ]trmes: 

// [ua] , 中 間 レ ベル 入力 の チャ ー ジ ・ ポ ンプ の 電流 値 , rp 用 $fdisp1ay(mcd, "ジッ タ 量 - - -RMS の 2 乗 値 : sgr, 
parameter CP=50.0: //[pF], LtPF で の キャ パシ タ 値 (フィ ル タ 定数 ) 
parameter C2=2.0: //[pF], チャ ー ジ ・ ポ ンプ 切り 替り 時 の (tmp 1/tmp_r2) ) : 

フィ ー ド スルー 緩和 用 キャ パシ タ 値 . rpg で 使用 $Edisp1ay (mcd, "VCO 生成 クロ ッ ク 半 周期 --- 最 小 値 : sq, 

parameter REP=30.0: //[Ohm], pF で の 抵抗 値 ( フ ィ ル タ 定数 ) 最大 値 : sd ,r3.tsh min gen c1k, 
parameter TS RK LPF=0.001: //[ns], ルン ゲ ク ッ タ 処理 の 周期 , ripg 3 .tsh max gen oc1k) : 
parameter N STP TNT LPF=0.064: //[ns], ルン ゲ ク ッ タ 計算 の $Fdisp1ay (mcd, "ルン ゲ ク ッ タ 計算 で 許容 誤差 を 満足 

ステ ッ プ 幅 の 初期 値 (rjpr 出力 の 更新 周期 の 初期 値 ) し な か っ た 回 数 : dr , エ 2.12.hit min h) 
parameter RELERR LPF=0.001: // 許 容 相対 エラ ー 項 (ルン ゲ ク ッ タ 計算) $Eo1oge (mcd) : 
parameter ABSERR LPF=0.00001, // 許 容 絶対 エラ ー 項 (ルン ゲ ク ッ タ 計算 ) $Eintgh: 
parameter CYCLE VCO=1, //VCO で の 出力 更新 時 間 (=TM_PRCS*CYCLE VCO) end 
parameter THETA TNTT=0.0, //[rag], VCO で の 初期 位相 


end 

// 測 定 回 路 (msjt) 用 パラ メー タ 

parameter CYCLE MJ=1 a1ways begin // 入 力 デー タ 生成 用 リフ ァ レ ンス ・ ク ロッ ク の 生成 
// 測 定 を 行う 際 の スキ ッ プ 数 . 1/Cycrg Mr が 測定 頻度 を 意味 する ref o1k=1'b0: #(tg reFo1k/2) , 

parameter T AVRG=300: ref co1k=1!D1: 后 (Eg reFo1k/2) : 
// [ns], strt msjt が エ 1- >H 後 の ジッ タ 平 均 値 を 求め る た め の 時 間 を 指す end 

parameter TT AVMJ=500: 
// [ns], ピー ク ・ ジ ッ タ や RMs ジ ッ タ 値 を 求め る た め の 時 間 を 指す a1ways begin 

parameter TS CLK=0.2: //VCOo チ ェ ー ン 回 路 制御 電圧 値 を 書き 込む た め の ク ロッ ク 
// [ns], リフ ァ レ ンス ・ ク ロッ ク の 周期 ms]t oc1k=1'b0, #5000: 

// <== 時 間 単 位 +00fs の 場合 に 変 

変 


生成 


defFparam T2 .TM PRCS=TM PRCS: mg 上 Cc1k=1!D1: 半 5000: 
defparam T2.K RTOT=K RTOT: // <== 時 間 単 位 100Es の 場合 に 
defparam 12.TD PD1=TD PD1: end 

defparam 12.TREF PD1=TRE PD1: 
defparam 12.K DIV PD1=K DTV PD1: a1way8 @ (ms]t c1k) begin 

defparam 12.TD FD=TD FD: //Vco チ ェ ー ン 回 路 入 力 (制御 ) 電圧 値 の 書き 込み 
defparam T2 .TRE FD=TREF EFD: $Fd1gp1ay (mcd, "を d gd",$ モ 1me, T2.T エ 3 .Yin 1): 
defparam 12.K DTV FD=K DTV EFD: //1 行 目 は 時 間 , 2 行 目 は vco 入力 

defparam T2.N STP TNT LPF=N STP TNT LPF: end 

defparam 12.RELERR LPF=RELERR LPE: 
defparam 12.ABSERR LPF=ABSERR LPE: endmodu1e 
defparam 12.TS RK LPF=TS RK LPF: 
defparam 12.CUR CP PD1=CUR CP PD1 
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